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We develop a new formulation of the functional renormalization group (RG) for interacting 
fermions. Our approach unifies the purely fermionic formulation based on the Grassmannian func- 
tional integral, which has been used in recent years by many authors, with the traditional Wilsonian 
RG approach to quantum systems pioneered by Hertz [Phys. Rev. B 14, 1165 (1976)], which at- 
tempts to describe the infrared behavior of the system in terms of an effective bosonic theory 
associated with the soft modes of the underlying fermionic problem. In our approach, we decouple 
the interaction by means of a suitable Hubbard-Stratonovich transformation (following the Hertz- 
approach), but do not eliminate the fermions; instead, we derive an exact hierarchy of RG flow 
equations for the irreducible vertices of the resulting coupled field theory involving both fermionic 
and bosonic fields. The freedom of choosing a momentum transfer cutoff for the bosonic soft modes 
in addition to the usual band cutoff for the fermions opens the possibility of new RG schemes. In 
particular, we show how the exact solution of the Tomonaga-Luttinger model (i.e., one-dimensional 
fermions with linear energy dispersion and interactions involving only small momentum transfers) 
emerges from the functional RG if one works with a momentum transfer cutoff. Then the Ward 
identities associated with the local particle conservation at each Fermi point are valid at every stage 
of the RG flow and provide a solution of an infinite hierarchy of flow equations for the irreducible 
vertices. The RG flow equation for the irreducible single-particle self-energy can then be closed and 
can be reduced to a linear integro-differential equation, the solution of which yields the result famil- 
iar from bosonization. We suggest new truncation schemes of the exact hierarchy of flow equations, 
which might be useful even outside the weak coupling regime. 



PACS numbers: 71.10.Pm., 71.10.Hf 



I. INTRODUCTION 



In condensed matter and statistical physics the renor- 
malization group (RG) in the form developed by Wilson 
and coauthorsi*2^ has been very successful. At the heart 
of this intuitively appealing formulation of the RG lies 
the concept of an effective action describing the physical 
properties of a system at a coarse grained scale. Con- 
ceptually, the derivation of this effective action is rather 
simple provided the theory can be formulated in terms 
of a functional integral: one simply integrates out the 
degrees of freedom describing short-wavelength fluctua- 
tions in a certain regime, and subsequently rescales the 
remaining degrees of freedom in order to compare the 
new effective action with the original one. Of course, in 
practice, the necessary functional integration can almost 
never be performed analytically, so that one has to resort 
to some approximate procedure^ However, if the elimi- 
nation of the short wavelength modes is performed in in- 
finitesimal steps, one can write down a formally exact RG 
flow equation describing the change of the effective action 
due to mode elimination and rescaling. The earliest ver- 
sion of such a functional RG equation has been derived 
by Wegner and Houghton. 5 Subsequently, many authors 
have derived alternative versions of the functional RG 
with the same physical content, using different types of 
generating functionals. In particular, the advantages of 



working with the generating functional of the one-particle 
irreducible vertices have been realized early on by Di Cas- 
tro, Jona-Lasinio, and Peliti^ and by Nicoll, Chang, and 
Stanleyi^ The focus of the above works was an accu- 
rate description of second-order phase transitions at finite 
temperatures, where quantum mechanics is irrelevant. 

In recent years, there has been much interest in quan- 
tum systems exhibiting phase transitions as a function 
of some nonthermal control parameter, such as pressure 
or density. In a pioneering paper, Hertz? showed how 
the powerful machinery of the Wilsonian RG can be gen- 
eralized to study quantum critical phenomena in Fermi 
systems. Technically, this is achieved with the help of 
so-called Hubbard-Stratonovich transformations, which 
replace the fermionic two-particle interaction by a suit- 
able bosonic field that couples to a quadratic form in the 
fcrmion operators^ The fermions can then be integrated 
out in a formally exact way, resulting in an effective ac- 
tion for the bosonic field. Of course, there are many pos- 
sible ways of decoupling fermionic two-body interactions 
by means of Hubbard-Stratonovich transformations. In 
the spirit of the usual Ginzburg-Landau- Wilson approach 
to classical critical phenomena, one tries to construct the 
effective bosonic theory such that the field can be iden- 
tified with the fluctuating order parameter, or its field 
conjugate. However, if the system supports additional 
soft modes that couple to the order parameter an 
attempt to construct an effective field theory in terms of 
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the order parameter alone leads in general to an effective 
action with singular and nonlocal vertices, so that the 
usual RG methods developed for classical phase tran- 
sitions cannot be applied. In this case, it is better to 
construct an effective action involving all soft modes ex- 
plicitly. However, for a given problem the nature of the 
soft modes is not known a priori, so that the explicit 
introduction of the corresponding degrees of freedom by 
means of a suitable Hubbard-Stratonovich transforma- 
tion is always based on some prejudice about the nature 
of the ground state and the low-lying excitations of the 
system. Recently, the breakdown of simple Ginzburg- 
Landau- Wilson theory has also been discussed in the con- 
text of quantum antiferromagnets by Senthil et alm& 

In order to have a completely unbiased RG approach 
for interacting fermions, one can also apply the Wilso- 
nian RG directly to the Grassmannian functional in- 
tegral representation of the partition function or the 
Green's functions of an interacting Fermi system^* In 
the past ten years, several groups have further devel- 
oped this method. 15 i 16 i 1 7i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 2fi i 2 7 In par- 
ticular, the mode elimination step of the RG transforma- 
tion has been elegantly cast into formally exact differen- 
tial equations for suitably defined generating functionals. 
These functional equations then translate to an infinite 
hierarchy of integrodifferential equations for the vertices. 
On a technical level, it is often advantageous to work 
with the generating functional of the one-particle irre- 
ducible vertices £L£ For classical field theories the exact 
RG flow equation for this generating functional has been 
obtained by Wetterich^ and by Morris,^ 9 For nonrel- 
ativistic fermions, the corresponding flow equation has 
been derived by Salmhofer and HonerkampAS and, in- 
dependently and in more explicit form, by Kopietz and 
BuscheiSS 

However, the purely fermionic formulation of the 
Wilsonian RG has several disadvantages. In order to ob- 
tain at least an approximate solution of the formally ex- 
act hierarchy of RG flow equations for the vertices, severe 
truncations have to be made, which can only be justified 
as long as the fermionic four-point vertex (i.e., the effec- 
tive interaction) remains small. Hence, in practice the 
fermionic functional RG is restricted to the weak cou- 
pling regime, so that possible strong coupling fixed points 
are not accessible within this method. Usually, one-loop 
truncated RG equations are iterated until at least one of 
the marginal interaction constants becomes large, which 
is then interpreted as a weak-coupling instability of the 
Fermi system in the corresponding channel^*A£*i& 

Within the framework of the Wilsonian functional 
RG, a consistent two-loop calculation has not been per- 
formed so far due to the immense technical difficulties 
involved. Such a calculation should also take into ac- 
count the frequency dependence of the effective inter- 
action and the self-energy corrections to the internal 
Green's functions. Note that in the work by Katanin 
and KampfS 6 . the frequency dependence of the effective 
interaction has been ignored, so that the effect of possible 



bosonic collective modes is not included in these calcu- 
lations. It is well known that two-loop calculations are 
more conveniently performed using the field-theoretical 
RG method, provided the physical problem of interest 
can be mapped onto a renormalizable field theory. Fer- 
raz and coauthors22i2i recently used the field-theoretical 
RG to calculate the single-particle Green function at the 
two loop level for a special two-dimensional Fermi sys- 
tem with a flat Fermi surface. Interestingly, they found 
a new non-Fermi liquid fixed point characterized by a 
finite renormalized effective interaction and a vanishing 
wave-function renormalization. The running interaction 
without wave-function renormalization diverges in this 
model, but the divergence is canceled by the vanishing 
wave-function renormalization such that the renormal- 
ized interaction remains finite. 

Another difficulty inherent in any RG approach to 
Fermi systems at finite densities arises from the fact 
that the true Fermi surface of the interacting many- 
body system is not known a priori. In fact, in dimen- 
sions D > 1 interactions can even change the symme- 
try of the Fermi surfacei 3 ^ 3 - In a perturbative approach, 
finding the renormalized Fermi surface is a delicate self- 
consistency problem 3 ^; if one starts from the wrong Fermi 
surface one usually encounters unphysical singularities^ 
So far, the self-consistent renormalization of the Fermi 
surface has not been included in the numerical analysis 
of the one-loop truncated fermionic functional RG. As 
shown in Refs. 

nana 

and l3ll the renormalized Fermi 
surface can be defined as a fixed point of the RG and can 
in principle be calculated self-consistently entirely within 
the RG framework. 

The progress in overcoming the above difficulties in- 
herent in the RG approach to fermions using a purely 
fermionic parametrization has been rather slow. In our 
opinion, this has a simple physical reason: the low ly- 
ing excitations (i.e., soft modes) of an interacting Fermi 
system consist not only of fermionic quasiparticles, but 
also of bosonic collective excitations. 36 The latter are 
rather difficult to describe with in the pu rely fermionic 
parametrizations used in Refs. ll5ll6lll^lTMl9ll20ll2lll22l 
I2^l2^l25ll2^l2^ . Naturally, a formulation of the functional 
RG where both fermionic and bosonic excitations are 
treated on equal footing should lead to a more convenient 
parametrization. Such a strategy seems also natural in 
light of the observation by Kirkpatrick, Belitz, and co- 
workers^ (see also Ref. lf2T) that an effective low-energy 
and long-wavelength action with well-behaved vertices is 
only obtained if the soft modes are not integrated out, 
but appear explicitly as quantum fields. 

Our aim in this work is to set up a functional renormal- 
ization group scheme that allows for a simultaneous treat- 
ment of fermionic as well as collective bosonic degrees of 
freedom. This is done by explicitly decoupling the inter- 
action via a Hubbard-Stratonovich transformation in the 
spirit of Hertz 9 and then considering the functional renor- 
malization group equations for the mixed field theory 
involving both fermionic and bosonic fields. This type 



3 



of approach has been suggested previously by Correia, 
Polonyi, and Ri chert }21 who studied the homogeneous 
electron gas by means of a gradient expansion of a func- 
tional version of a Callan-Symanzik equation. Here, we 
follow the more standard approach and derive a hierarchy 
of flow equations for the vertex functions of our coupled 
Fermi-Bose theory. A related approach has also been de- 
veloped by Wetterich and coauthors However, they 
discussed only a simple truncation of the exact hierarchy 
of RG flow equations involving an effective potential in 
a bosonic sector and a momentum-independent Yukawa 
coupling. They did not pay any attention to the problem 
of the compatibility of the RG flow with Ward identities, 
which play a crucial role for interacting fermions with 
dominant forward scattering. Our procedure bridges the 
gap between the purely bosonic approach to quantum 
critical phenomena by HertzS and the purely fermionic 
functional RG method developed in the past decade by 
several authorsi 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 Contrary 
to the approach by Hertz^ in our approach the fermionic 
degrees of freedom are still present, so that the feedback 
of the collective bosonic modes on the one-particle spec- 
tral properties of the fermions can be studied. 

The parametrization of the low-lying excitations of an 
interacting Fermi system in terms of collective bosonic 
fields is most natural in one spatial dimension, 40-41 where 
Fermi-liquid theory breaks down and is replaced by the 
Luttinger- liquid concept 4i An exactly solvable paradigm 
of a Luttinger liquid is the Tomonaga-Luttinger model 
(TLM)^LMJ1 consisting of fermions with exactly lin- 
ear energy dispersion and interactions involving only 
small momentum transfers. A description in terms of 
bosonic variables is well known to provide an exact so- 
lution for thermodynamic quantities as well as corre- 
lation functions^Siii^ One might then wonder if it is 
also possible to obtain the complete form of the corre- 
lation functions of the TLM entirely within the frame- 
work of the functional RG. An attempt^ to calculate 
the momentum- and frequency-dependent single-particle 
spectral function A(k,u>) of the TLM by means of an 
approximate iterative two-loop solution of the functional 
RG equations at weak coupling yields the correct behav- 
ior of A(±fcp-,w) known from bosonization (where kp is 
the Fermi momentum) , but yields incorrect threshold sin- 
gularities for momenta away from ±fcp. In this work we 
go considerably beyond Ref. |^] and show how the TLM 
can be solved exactly using our mixed fermionic-bosonic 
functional RG. A crucial point of our method is that the 
RG can be set up in such a way that the Ward identities 
underlying the exact solubility of the TL M 43 i 44 i 4 ^ 47 are 
preserved by the RG. This is not the case in the purely 
fermionic formulation of the functional RGi^ 7 . 

Very recently, Benfatto and Mastropietro^ 8 - also devel- 
oped an implementation of the RG for the TLM which 
takes the asymptotic Ward identities into account. How- 
ever, these authors did not introduce bosonic collective 
fields and did not attempt to calculate the exact single- 
particle Green's function. 



The rest of the paper is organized as follows. In Sec.lTTl 
we introduce the model, carry out the decoupling of the 
interaction and set up a compact notation which turns 
out to facilitate the bookkeeping in the derivation of the 
functional RG equations that is presented in detail in 
Sec. IIIII In Sec. IIVI we introduce a new RG scheme 
which uses the momentum transfer of the effective inter- 
action as a cutoff. We show that for linearized energy 
dispersion the resulting infinite hierarchy of flow equa- 
tions for the irreducible vertices involving two external 
fermion legs and an arbitrary number of external boson 
legs can be solved exactly by means of an infinite set of 
Ward identities. Using these identities, the flow equation 
for the irreducible self-energy can then be reduced to a 
closed linear integrodifferential equation, which can be 
solved exactly. In one dimension, we recover in Sec. lIVCl 
the exact solution of the Tomonaga-Luttinger model in 
the form familiar from functional bosonization 4 ^^. Fi- 
nally, in Sec. E] we summarize our results and give a brief 
outlook on possible further applications of our method. 
There are three appendices where we present some more 
technical details. In Appendix A, we use our compact 
notation introduced in Sec. ^ to discuss the structure 
of the tree expansion in our coupled Fermi-Bose theory. 
Appendix B contains a derivation of the skeleton dia- 
grams for the first few irreducible vertices of our theory 
using the Dyson-Schwinger equations of motion, which 
follow from the invariance of the functional integral with 
respect to infinitesimal shift transformations. Finally, in 
Appendix C we use the gauge invariance of the mixed 
Fermi-Bose action to derive a cascade of infinitely many 
Ward identities involving vertices with two fermion legs 
and an arbitrary number of boson legs. 



II. INTERACTING FERMIONS AS COUPLED 
FERMI-BOSE SYSTEMS 

In this section we discuss the Hubbard-Stratonovich 
transformation and set up a condensed notation to treat 
fermionic and bosonic fields on the same footing. This 
will allow us to keep track of the rather complicated di- 
agrammatic structure of the flow equations associated 
with our coupled Fermi-Bose system in a very efficient 
way. A similar notation has been used previously in 
Refs.dandli! 



A. Hubbard-Stratonovich transformation 

We consider a normal fermionic many-body system 
with two-particle density-density interactions. In the 
usual Grassmannian functional integral approach 1 ^ the 
grand-canonical partition function and all (imaginary)- 
time-ordered Green's functions can be represented as 
functional averages involving the following Euclidean ac- 
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tion: 



Siat[i>,ip] 



S Q [>P,i>}+S int [^,i}}, (2.1) 

= V\ I $K a [-*U + buv]ll>Ko , (2-2) 
2 J R PKaPRo' ' 



(2.3) 



where the composite index A" = (iw,k) contains a 
fermionic Matsubara frequency as well as an ordinary 
wave vector k. Here the energy dispersion £ker = £ko- — M 
is measured relative to the chemical potential n, and 
are some momentum-dependent interaction parameters. 
Throughout this paper, labels with an overbar refer to 
bosonic frequencies and momenta, while labels without 
an overbar refer to fermionic ones. We have normalized 
the Grassmann fields ipKv and tpKa such that the inte- 
gration measure in Eq. (|2.1(l is 



1 



x - /3,V->oo 



K 



U.K. 



duj d D k 
2^ {2tt) d 



(2.4) 



where (3 is the inverse temperature and V is the volume 
of the system. The Fourier components of the density 
are represented by the following composite field: 



PRa 



K 



K+K.a ? 



(2.5) 



which implies p Ka — p_R a . The discrete index a is for- 
mally written as a spin projection, but will later on also 
serve to distinguish right and left moving fields in the 
Tomonaga-Luttinger model. This is why a dependence 
of the dispersion on a has been kept. 

The interaction is bilinear in the densities and can 
be decoupled by means of a Hubbard-Stratonovich 
transformation^ The interaction is then mediated by 
a real field tp and the resulting action reads as 

S$,il>M = S $,il>]+S [<p] + SS,il>M, (2-6) 
where the free bosonic part is given by 



sow - It, /[/t 1 ] 

era' 

and the coupling between Fermi and Bose fields is 



K 



{ T / ^K+K^KafRa- ( 2 - 

~ JK JK 



The Fourier components of a real field satisfy <p* Ka — 
tp-Ko-- For the manipulations in the next section it will 
prove advantageous to further condense the notation and 
collect the fields in a vector $ = (tp, ^j, ip). The quadratic 



part of the action can then be written in the symmetric 
form 



s m = sS, i>] + SoM = -\ (*, [Go]- 1 $ 



$ Q [G ];i So? , 



(2.9) 



where Go is now a matrix in frequency, momentum, spin, 
and field-type indices, and a is a "super label" for all of 
these indices. The symbol J denotes integration over 
the continuous components and summation over the dis- 
crete components of a. The matrix Gq 1 has the block 
structure 



Go 1 



o cfoT 





1 




(2.10) 







whereS . £ = — 1 and Go and Fq are infinite matrices 
in frequency, momentum, and spin space, with matrix 
elements 

[Go\Ka,K'a' = $K, K> 8 a a' G ,er ( K) , (2.11) 

[Fo]k*,k'*> = t K +K>,oFo,™>(K) , (2.12) 



where 



Fo,crcr> (K) = Si 



a a 
k 



(2.13) 
(2.14) 



The Kronecker Sk.k' — PV&u>,uj'&\t,\t' appearing in 
Eqs. I|2.11I2.12|I is normalized such that it reduces to 
Dirac 8 functions 8k, K' ~ * (2tt) d+1 8(u — u>' )«5( D )(k-k') 
in the limit (3, V — > oo. Note that the bare interaction 
plays the role of a free bosonic Green's function. For later 
reference, we note that the inverse of Eq. (|2.1U|) is 



G 
Go = ( CGJ 
-F 



(2.15) 



and that the transpose of Go satisfies 

G^ = ZG = G Z , (2.16) 

where the "statistics matrix" Z is defined by 

[Zjcra' = <5 QQ 'Ca ■ (2-17) 

Here, £<* = — 1 if the supcrindex a refers to a Fermi field, 
and Ca = 1 if a labels a Bose field. 

B. Generating functionals 

1. Generating functional of connected Green's functions 

We now introduce sources J a and define the generating 
functional Q[J] of the Green's functions as follows: 

G[J] = e e °W = _L [ m e -s -Sx+(./,*) . ( 2 .18) 

Zq J 
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Here, Q C [J] is the generating functional for connected 
Green functions and the partition function Z of the non- 
interacting system can be written as the Gaussian inte- 
gral 



Z a = JD<I> e~ So 
Let us use the compact notation 



(J,$)= f J Q $ C 

J a 



(2.19) 



(2.20) 



Conventionally, the source terms for fields of different 
types are written out explicitly in the formic 

Y] / jKalpKa + Y" / 4>K*jl<a + Y" / Jr^Kct ■ 

(2.21) 



K 



A comparison between Eq. (|2.2U|) and Eq. 12.21fl shows 
that the sources in the compact notation are related to 
the standard ones by J = (J, (j, J* ) . The connected n- 

line Green's functions Gc"ai...a„ are then defined via the 
functional Taylor expansion 

Gc[J] = £ ~ } I ■■ I ^ c %... a J ai ••• J an , (2.22) 
„=o n - Ja i Ja « 



implying 



SJ a 



, SJ ai 



(2.23) 



j=o 



In particular, the exact Green's function of our interact- 
ing system is given by 



SJ a SJ a i 



T (2) 



(2.24) 



j=o 



which we shall write in compact matrix notation as 



G 



5JSJ 



j=o 



G 
CG T 
—F 



(2.25) 



For the last equality, it has been assumed that no sym- 
metry breaking occurs. Thus, G has the same block 
structure as the noninteracting Go in Eq. I|2.15|l so that, 
similarly to Eq. (|2.1ti|) . 



ZG = GZ 



(2.26) 



In the noninteracting limit (Si 0) one easily verifies by 
elementary Gaussian integration that the matrix G given 
in Eq. I|2.24|l reduces to G , as defined in Eq. I|2.15|l . The 
self-energy matrix also has the same block structure as 



the inverse free propagator. Dyson's equation then reads 
as 



G- 1 = Go 1 



(2.27) 



where the matrix S contains the one-fermion-line irre- 
ducible self-energy T, a (K) and the one-interaction- line 
irreducible polarization II CT (^") in the following blocks: 

o cp] T o \ 

£ , (2.28) 

o o n / 

where 

[£W,k>*- = 5 K>K '5 aal ^ a {K') , (2.29) 
= 5k + k',o^^{K'). (2.30) 

These matrices are spin-diagonal because the bare cou- 
pling Si[ip,ip,(p] between Fermi and Bose fields in 
Eq. H2.8fl is diagonal in the spin index. The blocks of 
the full Green's function matrix G in Eq. I|2.25|l contain 
the exact single-particle Green's function and the effec- 
tive (screened) interaction, 



1 F ]k<j,K'<t> 



5K,K'Saa'G a (K) , 

S K+R',n F aa -(K) , 



with 



G„{K) = [G^m-^iK)]- 1 
-l 



F aa ,(K) = 



n 



Ka,-Ka' 



(2.31) 
(2.32) 



(2.33) 
(2.34) 



2. Generating functional of one-line irreducible vertices 

Below, we shall derive exact functional RG equations 
for the one-line irreducible vertices of our coupled Fermi- 
Bose theory. The diagrammatic perturbation theory con- 
sists of both fermion and boson lines. We require irre- 
ducibility with respect to both types of lines and call this 
one-line irreducibility. One should keep in mind that 
a boson line represents the two-body electron-electron 
interaction which is screened by zero-sound bubbles for 
small momentum transfers. This means that in fermionic 
language our vertices are not only one-particle irreducible 
but are also approximately two-particle irreducible in the 
zero-sound channel in the sense that particle-hole bubbles 
are eliminated in favor of the effective bosonic propaga- 
tor. In order to obtain the generating functional of the 
corresponding irreducible vertices, we perform a Legen- 
dre transformation with respect to all field components, 
introducing the classical field^i 



SGc 

SJ a 



(2.35) 
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After inverting this relation for J = J[$] we may calcu- 
late the Legendre effective action 



£[$] = (J[$],$) -&[■/[<&]] 
From this we obtain 

5C 



(2.36) 



J a — Co: 



(2.37) 



which we may write in compact matrix notation as 

.sc 



J 



(2.38) 



In this notation the chain rule simply reads as 



_S_ 

6$ 



8W£„ S 
- Zi - 



5<P5<P 5J 

Applying this to both sides of Eq. (|2.35|) we obtain 



(2.39) 



_ 6$ 
1 ~ M> 



5<i>5<f> 8J8J 
For vanishing fields $ and J this yields 

6^C 



<5<M<I> 



= ZG 1 = fG 



liT 



(2.40) 



(2.41) 



The advantage of our compact notation is now obvious: 
the minus signs associated with the Grassmann fields 
can be neatly collected in the "statistics matrix" Z. If 
the Grassmann sources are introduced in the conven- 
tional way^fi the minus signs generated by commuting 
two Grassmann fields are distributed in a more compli- 
cated manner in the matrices of second derivatives^^ 
From Eq. I|2.41|l it is evident that we need to subtract 
the free action from £[$] to obtain the generating func- 
tional for the irreducible vertex functions, 

= £[*] - Soffl = C[<f} + [G^ . (2.42) 
Then we have, using the Dyson equation (|2.27|) . 



<5$<5$ 



8^C 



<s>=o 



5<$>5<$> 
-\G 



-liT 



$=0 



IG 



o J 

-liT 
J 



(2.43) 



In general, the one-line irreducible vertices are defined as 
coefficients in an expansion of r[$] with respect to the 
fields, 



n=0 U - Jq i 



p(n) 



(2.44) 



The vertices T^ n ' have the same symmetry with respect 
to interchange of the indices as the monomial in the fields, 
i.e., the interchange of two neighboring Fermi fields yields 
a minus sign. Graphically, we represent the vertices 1^™) 



n-1 




FIG. 1: Graphical representation of the symmetrized one-line 
irreducible n-point vertex defined via Eq. 12.441 . Because for 
fermions the order of the indices is important (exchanging two 
neighboring fermion legs will generate a minus sign) the circles 
representing the irreducible vertices have an arrow that points 
to the leg corresponding to the first index; subsequent indices 
are arranged in the order indicated by the arrow. External 
legs denote either outgoing fermions (tp), incoming fermions 
(ip), or bosons (ip). 



by an oriented circle with n external legs, as shown in 
Fig.Hl With the definition $F%^ and Eq. fFZty we have 

T^', — [S] a r a . The fact that also the higher-order ver- 
tices r(") defined in Eq. (|2.44|) are indeed one- line ir- 
reducible (i.e., cannot be separated into two parts by 
cutting a single fermion line or a single interaction line) 
can be shown iteratively by generating a tree expansion 
from higher-order derivatives of Eq. (|2.4U|I . We show this 
explicitly in Appendix A. 



III. FUNCTIONAL RG FLOW EQUATIONS 
FOR ONE-LINE IRREDUCIBLE VERTICES 

In this section we derive exact RG flow equations for 
the generating functional of the one-line irreducible ver- 
tices of our coupled Fermi-Bose theory. We also classify 
the various vertices of the theory based on their scal- 
ing dimensions and propose a new truncation scheme in- 
volving the building blocks of the skeleton diagrams for 
fermionic and bosonic two-point functions. 



A. Cutoff schemes 

Since the interaction now appears as a propagator of 
the field (p, it is possible to introduce a momentum- 
transfer cutoff in the interaction on the same footing as 
a bandwidth cutoff. A bandwidth cutoff restricts the rel- 
evant fermionic degrees of freedom to the vicinity of the 
Fermi surface, and appears to be most natural in the 
RG approach to fermions in one spatial dimension^ In 
higher dimensions, the Wilsonian idea of eliminating the 
degrees of freedom in the vicinity of the Fermi surface 
is implemented by defining for each momentum k an as- 
sociated hp by means of a suitable projection onto the 
Fermi surface^ and then integrating over fields with mo- 
menta in the energy shell vqA < |ek — £k F | < vqAq, where 
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Ek is the energy dispersion in the absence of interactions. 
Here is some suitably defined velocity (for example, 
some average Fermi velocity) , which we introduce to give 
A units of momentum. We shall refer to voA as band- 
width cutoff. Formally, we introduce such a cutoff into 
our theory by substituting for the free fermionic Green's 
function in Eq. l|2.I3f) 



GoA K ) 



Q(A<D K <A )G 0ja (K) 
0(A < D K < A ) 



•6 



with 



D K = |ek - £k F |/wo 



(3.1) 



(3.2) 



Here 0(A < x < Ao) = 1 if the logical expression in the 
brackets is true, and 0(A < x < A ) = if the logical 
expression is false. Ambiguities associated with the sharp 
O-function cutoff can be avoided by smoothing out the 
functions and taking the sharp cutoff limit at the end 
of the calculation^ In order to construct a consistent 
scaling theory, the kp in Eq. (|3.2() should refer to the 
true Fermi surface of the interacting system, which can in 
principle be obtained self-consistently from the condition 
that the RG flows into a fixed pointi 20 i 22 

The above bandwidth-cutoff procedure has several dis- 
advantages. On the one hand, for any finite value of the 
cutoff parameter vqA the Ward identities are violated. 27 
Moreover, the RG flow of two-particle response functions 
probing the response at small momentum transfers (such 
as the compressibility or the uniform magnetic suscep- 
tibility) is artificially suppressed by the bandwidth cut- 
off. To cure the latter problem, various other parame- 
ters have been proposed to serve as a cutoff for the RG, 
such as the temperatureiS. or even the strength of the 
interaction^ While for practical calculations these new 
cutoff schemes may have their advantages, the intuitively 
appealing RG picture that the coarse grained parameters 
of the renormalized theory contain the effect of the de- 
grees of freedom at shorter length scales and higher en- 
ergies gets somewhat blurred (if not completely lost) by 
these new schemes. 

The above mentioned problems can be elegantly 
avoided in our mixed Fermi-Bose theory if we work with 
a momentum cutoff in the bosonic sector of our theory, 
which amounts to replacing in Eq. I|2.14[l . 



Ft 



O.crcr 



(K) 



Q(A<D R <A )Fo !(T(T ,(K) 

= 0(A<D R < A )f™' , (3.3) 



where 



(3.4) 



Keeping in mind that the bosonic field mediates the ef- 
fective interaction, it is clear that A is a cutoff for the 
momentum transfer of the interaction. This is precisely 
the same cutoff scheme employed in the seminal work 



by Hertz ^ who discussed also more general frequency- 
dependent cutoffs for the labels of the bosonic Hubbard- 
Stratonovich fields, corresponding to more complicated 
functions D R than the one given in Eq. (|3.4|) . Moreover, 
in the exact solution of the one-dimensional Tomonaga- 
Luttinger model (abbreviated here as TLM, as already 
defined above) by means of a careful application of the 
bosonization method^ the maximal momentum trans- 
fered by the interaction appears as the natural cutoff 
scale. 

In our RG approach we have the freedom of choos- 
ing both the bandwidth cutoff vqA and the momentum- 
transfer cutoff A independently. In particular, we may 
even choose to get rid of the bandwidth cutoff completely 
and work with a momentum-transfer cutoff only. In this 
work, we shall show that if the interaction involves only 
small momentum transfers, then the pure momentum- 
transfer cutoff scheme indeed regularizes all infrared sin- 
gularities in one dimension. Moreover and most impor- 
tantly, introducing a cutoff only in the momentum trans- 
fer leads to exact RG flow equations that do not violate 
the Ward identities responsible for the exact solubility of 
the TLM. Given this fact, it is not surprising that we can 
solve the infinite hierarchy of RG flow equations exactly 
and obtain the exact single-particle Green's function of 
the TLM within the framework of the functional RG. 



Flow equations for completely symmetrized 
vertices 



With the substitutions (|3.1|l and (|3.3|) the noninteract- 
ing Green's function Go, and hence all generating func- 
tionals, depend on the cutoff parameter A. We can now 
follow the evolution of the generating functionals as we 
change the cutoff. The differentiation of Eq. (|2.18l) with 
respect to A yields for the generating functional of the 
Green's functions, 



d A G 



For the connected version Q C [J] — \nQ[J], we obtain 



SJSJ 



d A In Zo. (3.6) 



In the derivation of flow equations for £ or T [see 
Eqs. i|2.36[l and (|2.42|) ]. we should keep in mind that in 
these functionals the fields <E> are held constant rather 
than the sources J. Hence, Eq. (|2.36(1 implies 



d A c[<f} = -d A g c [j}\ J= 



j Am 



(3.7) 
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FIG. 2: Graphical representation of Eq. 113. 1611 for n — 2, 
describing the flow of the totally symmetric two-point ver- 
tex. Empty circles with G and G denote the exact matrix 
propagator G and the single-scale propagator G defined in 
Eqs. 12.24H and H3.10H . respectively. 



Using this and Eq. I|3.6fl we obtain for the functional 

r[#] =£[*]- s m, 



9aT= - -Tr ^[Go 1 ] 



SJSJ 



d A In Z 



(3.8) 



To derive a closed equation for T, we express the matrix 

<5 (2) Sc 



SJSJ in terms of derivatives of T using Eq. (|A4I) . After 
some rearrangements we obtain the exact flow equation 
for the generating functional of the one-line irre- 

ducible vertices, 



— Tr 
2 



ZG T U T {l-G T U T }~ 



ZG^S T {1-G^E T } , (3.9) 



where the matrix U[<£>] is the field-dependent part of 
the second functional derivative of r[$], as defined in 
Eq. ijXTJ. For convenience wc have introduced the single- 
scale propagator G as 

G = — G9a[G " 1 ]G = [1-GoS]- 1 (d A G ) [l-SCop 1 , 

(3.10) 

which reduces to Go = <9aGo in the absence of interac- 
tions. The matrix G has the same block structure as the 



matrix G in Eq. I|2.25|> . We denote the corresponding 

blocks by G and F. In the limit of a sharp 0-function 
cutoff 2 - the blocks of the single scale propagator are ex- 
plicitly given by 

[G] K *,K>*> = &K,K'&vv>G a (K) , (3.11) 

[hue,** = tii+K>,oF^(K), (3.12) 



with 

G a {K) 



S(A - D K ) 



iu) - (,] 



-S(A-D R ) 



n 



Ko-Ko' 



(3.13) 
,(3.14) 



where on the right-hand side of Eq. (|3.14|) it is understood 
that the 8-function cutoff should be omitted from the 
matrix elements of F§. 

The second line in Eq. (|3.9(l does not depend on the 
fields any longer and therefore represents the flow of the 
interaction correction r^ ) to the grand-canonical poten- 
tial, 



9 A r(°) 



i. 



-Tr 



ZG^l-G^y 1 



(3.15) 



Since we have already dropped constant parts of the ac- 
tion in the Hubbard-Stratonovich transformation, we will 
not keep track of in the following. 

The first line on the right-hand side of Eq. I|3.9|l gives 
the flow of one-line irreducible vertices. We can generate 
a hierarchy of flow equations for the vertices by expand- 
ing both sides in powers of the fields. On the left-hand 
side, we simply insert the functional Taylor expansion 
(I2.44|) of r[$], while on the right-hand side we substitute 
the expansion of U[$] given in Eq. (|A7|) . For a compar- 
ison of the coefficients on both sides, the right-hand side 
has to be symmetrized with respect to external lines on 
different vertices. We can write down the resulting infi- 
nite system of flow equations for the one-line irreducible 
vertices r'"' with n > 1 in the following closed form: 



1 ' v v 

^A^ai,. ..,a„ = ""o ^ / ^ ] ^n,mi+...+mi ^ai,...,a mi ;a mi + 1 ,...,a mi + m2 ;...-,a mi + ... + m[1 + 1 ,...,a n | 
l — l mi,...,m ; - 1 

xTr [zfiMZtl'G'r^^^G 2, ■ ■•rt' 1 + + 2) L,_ 1+1 ,...^] } ■ ( 3 - 16 ) 

I 



Here the matrices T^^.J m are given in Eq. i|A8() and the 
symmetrization operator S is defined in Eq. (|A11(I . The 
effect of S is rather simple: it acts on an expression al- 
ready symmetric in the index groups separated by semi- 
colons to generate an expression symmetric also with re- 
spect to the exchange of indices between different groups. 



From one summand in Eq. i|3.16[) the symmetrization op- 
erator S thus creates n!/(mi! ■ . . . • to/!) terms. 

Figures. [21 and [3] show a graphical representation of 
the flow of the vertices and T^ 3 ). With the graphi- 
cal notation for the totally symmetric vertices introduced 
in Fig. ^ all the signs and combinatorics have a graphi- 



9 




FIG. 3: Graphical representation of Eq. 13.161 for n = 3, de- 
scribing the flow of the totally symmetric three-point vertex. 

cal representation. In the next section we will leave the 
shorthand notation and go back to more physical ver- 
tices, explicitly exhibiting the different types of fields. All 
this can be done on a graphical level and involves only 
straightforward combinatorics. In this sense the deriva- 
tion of higher flow equations is at the same level of com- 
plexity as ordinary Feynman graph expansions. 

C. Flow equations for physical vertices 

Usually, the generating functional T[tp, tp, tp] is ex- 
panded in terms of correlation functions that are not 



symmetrized with respect to the exchange of legs involv- 




FIG. 4: Pictorial dictionary to translate graphs involving to- 
tally symmetrized vertices to ones involving physical vertices, 
which are only symmetrized within helds of the same type. 
The diagrams on the right-hand sides of the last four lines 
represent G, G, F and F respectively. 

ing different types of fields. If we explicitly display mo- 
mentum and frequency conservation, such an expansion 
reads as 



T[ip,ip,tp] 



oo oo 

EE 

n— m— 



(n!) 2 m! 



K' n a' n JK 1 a 1 



K m a„ 



^K' 1 + ...+K' n ,K 1 + ...+K n +K 1 - 



K„ 



XT 



(2n,m) 



(K 1 a 1 , . . . , K' n a' n ; K\<j\, . . . , K n a n ; K\a\, . . . , K m a m ) 



' ^K' n <j' n i>K iai ■ ■■■■ 1pK n a n 'fiK 1 a 1 



(3.17) 



Diagrammatically, we represent a physical vertex r( 2 "' m ) 
involving 2n external fermion legs and m external boson 
legs by a triangle to emphasize that our theory contains 
three types of fields, see Fig. 0] We represent a leg asso- 
ciated with a tjj field by an arrow pointing outward, a leg 
for i\) by an arrow pointing inward, and a leg for <p with 
a wiggly line without an arrow. Recall that our Bose 
field is real because it couples to the density, so that it 
should be represented graphically by an undirected line. 
On the contrary, for a propagator G or G, the field ip 
is represented by an arrow pointing inward and ip by an 
arrow pointing outward. Apart from the energy- and 



momentum-conserving delta function, the totally sym- 
metric vertices defined by the expansion (|2.44[) coincide 
with the nonsymmetric ones in Eq. H3.17(l for the same or- 
der of the indices. We can therefore obtain the flow equa- 
tions for the nonsymmetric vertices by choosing a definite 
realization of the external legs and by carrying out the 
intermediate sums over the different field species, i.e., by 
drawing all possible lines in the intermediate loop (two 
possible orientations of solid lines or one wiggly line). On 
the right-hand side one then has to appropriately order 
all the legs on the vertices, keeping track of signs for the 
interchange of two neighboring fermion legs. Having done 
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FIG. 5: Flow of the irreducible fermionic self-energy. The 
diagrams are obtained from the diagrams shown in Fig. El 
by specifying the external legs to be one outgoing and one 
incoming fermion leg. 




FIG. 6: Flow of the irreducible polarization, obtained from 
the totally symmetric diagram in Fig. El by setting both ex- 
ternal legs equal to boson legs. Note that each closed fermion 
loop gives rise to an additional factor of £ = — 1. 



so, we can use the pictorial dictionary in Fig.EJto obtain 
diagrams involving the physical correlation functions. In 
this way, we obtain from the diagram for the completely 
symmetric two-point vertex shown in Fig- Othe diagram 
for the fermionic self-energy in Fig. El as well as the di- 
agram for the irreducible polarization shown in Fig. El 
Moreover, if we specify the external legs in the diagram 
for the completely symmetric three-legged vertex shown 
in Fig.EJto be two fermion legs and one boson leg, we ob- 
tain the flow equation for the three-legged vertex shown 
in Fig. El 

The flow equation for the vertex correction in Fig. [7| 
looks very complicated, so that at this point the reader 
might wonder how in one dimension we will be able to 
obtain the exact solution of the TLM using our approach. 
We shall explain this in detail in Sec. lIVI but let us antic- 
ipate here the crucial step: obviously all diagrams shown 
in Figs. El El and [7| can be subdivided into two classes: 
those involving a fermionic single-scale propagator (the 
slash appears on an internal fermion line), and those with 
a bosonic single-scale propagator (with a slash on an in- 
ternal boson line). At this point we have not specified the 




FIG. 7: Flow of three-legged vertex with two fermion legs and 
one boson leg, obtained as a special case of the diagram in 
Fig. 

cutoff procedure, but as mentioned in the Introduction, 
in Sec. IIVI we shall work with a bosonic cutoff only. In 
this case all diagrams with a slash attached to a fermion 
line should simply be omitted. This is the crucial simpli- 
fication which will allow us to solve the hierarchy of flow 
equations exactly. Let us proceed in this section without 
specifying a particular cutoff procedure. 



D. Rescaling and classification of vertices 

In order to assign scaling dimensions to the vertices, we 
have to define how we rescale momenta, frequencies, and 
fields under the RG transformation. The rescaling is not 
unique but depends on the nature of the fixed point we 
are looking for. Let us be general here and assume that in 
the bosonic sector the relation between momentum and 
frequency is characterized by a bosonic dynamic expo- 
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nent (this is the exponent z introduced by Herta 9 -), 
while in the fermionic sector the corresponding dynamic 
exponent is z$. Rescaled dimensionless bosonic momenta 
q and frequencies e are then introduced as usual 9 - 

q = k/A , e = Lj/n A , oc A 2 ^ . (3.18) 

For convenience, we choose the factor Cl A such that it has 
units of energy; e is then dimensionless. 

The proper rescaling of the fermionic momenta is not 
so obvious. Certainly, all momenta should be measured 
with respect to suitable points kp on the Fermi surface. 
One possibility is to rescale only the component ku = 
(k — kp) • vj? of a given momentum that is parallel to the 
local Fermi surface velocity vf (and hence perpendicular 
to the Fermi surface) ji^S where Vp is a unit vector in 
the direction of vjr. Unfortunately, in dimensions D > 1 
this leads to rather complicated geometric constructions, 
because in a fixed reference frame the component kn to be 
rescaled varies for different points on the Fermi surface. 
However, if the initial momentum-transfer cutoff Ao in 
Eq. (|3.3(l is small compared with the typical radius of the 
Fermi surface, the initial and final momenta associated 
with a scattering process lie both on nearby points on 
the Fermi surface. In this case it seems natural to pick 
one fixed reference point kp^ on the Fermi surface, and 
then measure all fermionic momentum labels k^ and k^ in 
r (2«,m) (K / ! ... iK ' n ;K 1 ,...,K n ;K 1 ,...,K m ) relative to 
this kp.o.. Here, the index a labels the different points on 
the Fermi surface, for example, in one dimension a = ±1, 
with kp,±\ = ikp. We then define rescaled fermionic 
momenta q and frequencies e as follows: 

q= (k-kp^VA , e = uj/n A , Q\ oc A 2 * . (3.19) 

The factor Q\ should again have units of energy such 
that e is dimensionless. Iterating the usual RG steps 
consisting of mode elimination and rescaling, we then 
coarse grain the degrees of freedom in a sphere around 
the chosen point ki? CT . Because by assumption the max- 
imal momentum transfer mediated by the interaction is 
small compared with |kj? i<T |, the fermionic momenta ap- 
pearing in T^ m \K' x , . . '. , K' n , K u . . . , K n ; K u . . . , K m ) 
are all in the vicinity of the chosen kp\ CT . This property is 
also responsible for the approximate validity of the closed 
loop theorem for interacting fermions with dominant for- 
ward scattering in arbitrary dimensions. 4 ^' 4 - 7 ' 49 

Apart from the rescaling of momenta and frequencies, 
we have to specify the rescaling of the fields. As usual, we 
require that the Gaussian part S$[<fr] = So{4>, ij)] + So[<p] of 
our effective action is invariant under rescaling. For the 
fermionic part this is achieved by defining renormalized 
fields ipQa in D dimensions via 

/ Z \ V2 „ 

where Z is the fermionic wave-function renormaliza- 
tion factor and Q — (q, ie) denotes the rescaled 



fermionic momenta and Matsubara frequencies as de- 
fined in Eq. I|3.19|l . With this rescaling the wave- function 
renormalization and the Fermi velocity have a vanishing 
scaling dimension (corresponding to marginal couplings), 
while the momentum- and frequency-independent part of 
the self-energy is relevant with scaling dimension + 1 ; see 
Ref. |23- Analogously, we find that the bosonic Gaussian 
part of the action is invariant under rescaling if we ex- 
press it in terms of the renormalized bosonic field <pQ a 
defined by 

( z \ 1/2 
^=(a^J (3 ' 21) 

where Z is the bosonic wave-function renormalization 
factor, Q — (q, ie) denotes the rescaled bosonic momenta 
and Matsubara frequencies defined in Eq. I|3.18|l . and vq 
denotes the noninteracting density of states at the Fermi 
surface. We introduce the factor of v$ for convenience 
to make all rescaled vertices dimensionless. By construc- 
tion Eq. (|3.21l) assigns vanishing scaling dimensions to 
the bare interaction parameters f^ 7 , corresponding to 
marginal Landau interaction parameters. 

Expressing each term in the expansion of the generat- 
ing functional T[ip,ip,(p] given in Eq. (|3.17l) in terms of 
the rescaled variables defined above and using the fact 
that r is dimensionless, we obtain the scaling form of the 
vertices. Omitting for simplicity the degeneracy labels tr, 
and assuming z$ < z v ^ we define the rescaled vertices, 

f ; ^(Qlj ■ • ■ j Qni Qli ■ • • i Qn] Qli ■ ■ ■ j Qrn) = 

u -m/2 AD{n - 1+m/ 2) il -i n rn/2 z n 2m /2 

xT^ n ' m) (K[, K' n ; Ki, . . . . A,,: A'; K m ) , (3.22) 

where we have to exclude the cases of purely bosonic 
vertices {n = 0) as well as the fermionic two-point vertex 
(i.e., the rescaled irreducible self-energy, corresponding to 
n = 1 and m = 0), which both need separate definitions. 
For the purely bosonic vertices (n = 0) we set 

tf' m \Q u ...,Q m ) = ^ m/2 (A D n A )- 1+m / 2 z^ 2 

x Yf m \K x ,...,K m ) , (3.23) 

while for the fermionic two-point vertex we should sub- 
tract the exact fixed point self-energy S*(kir !(r , iO) at the 
Fermi-surface reference-point ki? CT and for vanishing fre- 
quency as a counterterm^2i2S 

if < 0) (Q;Q) ee E,(Q) = A p(jf) _ *0)] . 

' 'a 

(3.24) 

If necessary, the counterterm S»(ki? j(T , iO) can be re- 
constructed from the condition that the constant part 
f; = £j(0) of the self-energy flows into an RG fixed 
point We consider the rescaled vertices to be func- 
tions of the logarithmic flow parameter I = — ln(A/Ao). 
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Introducing the flowing anomalous dimensions associated we can then write down the flow equations for the 
with the fermionic and bosonic fields, rescaled vertices. Omitting the arguments, we obtain for 

n > 1 the flow equation- 5 

m = -dilnZ , fj t = -dilnZ, (3.25) 



diT\ 



1 1 r\ t~\ lib r\ 

ni , . m V— o „ o „ \ - o 

( I " )D + : mill - -(£> + ^) - n»» - — '/; - • t^t • Qi ■ ^) Qi • ^ 



2 v " 2 

p(2n,m) 
+ 1 I ' 

where z m ; n = minj^, z^p}. For n = 0we obtain from Eq. (|3.23|l . 



i(2n,m) 



d[ f(0,m) = 



III Q 



=i(0,' 



+ r 



(0,m) 



(3.26) 



(3.27) 



where we have introduced the notation 

Q — = q. V q + z^ e— , (3.28) 

s *- v * + *4- (3 - 29) 

The inhomogeneities in Eqs. H3.26|) and (|3.27|) are given 
by the rescaled version of the right-hand sides of the flow 
equations for the unrescaled vertices, i.e., for n > 1, and 

Zip ^ Zcp^^ 

' \Qli ■ ■ ■ ! Qn'i Qlt ■ ■ ■ i Qnl Ql; ■ ■ ■ > Qrn) = 
!y -m/2 A r.(n-l+m/2) n -lQm/2 2 „^ ro /2 

x[-Ad A ri 2n < m) (K';^;^})] , (3.30) 

and for n = 0, 

ri°' m) (Qi,...,g ro ) = ^ m/2 (A D n A )- 1+m / 2 z^ 2 

x[-A9 A rf m) (X!,...,X m )]. (3.31) 

By properly counting all factors it is then not difficult to 
see that the explicit expressions for the inhomogeneities 
in Eqs. (|3.3U|I and l|3.31|l can be simply obtained from 
their unrescaled counterparts by replacing all vertices 
and propagators with their rescaled analogs, where the 
rescaled propagators are defined by 

G(K) = -£-G(Q) , F{K) = -F{Q), (3.32) 

and the corresponding rescaled single-scale propagators 
arc defined via 

A6(K)=~d(Q) ) AF(K) = --F(Q). (3.33) 

11 A V Q 

From Eqs. I|3.26[l and (|3.27|) we can read off the scal- 
ing dimensions of the vertices: the scaling dimension of 



p(2n,m) m £) dimensions is 

£j(2n,m) _ 

f (1 - ri)D + z min - (D + z v )m/2 for n > 1 
\ (D + ^)(l-m/2) forn = ' 

(3.34) 

In the particular case of the Tomonaga-Luttinger model, 
where D = 1 and z$ = z v = 1, we have £)( 2rl > m ) = 
2-n-m. Hence, in this case f( 2 <°)(Q = 0) and f (0 >^ are 
relevant with scaling dimension +1, while T^ 4 '°\Qi = 0) 
and r' 2,1 5(Qj = Qi = 0) are marginal. All other ver- 
tices are irrelevant. Of course, the linear terms in the 
expansion of f( 2 '°)(Q;Q) for small Q are also marginal. 
These terms determine the wave-function renormaliza- 
tion factor Z and the Fermi velocity renormalization vi, 
see Eqs. (|4.2U|I and (|4.25|) below. Note that for short- 
range interactions the dispersion of the zero-sound mode 
is linear in any dimension, 47 Hence, as long as the density 
response is dominated by the zero sound mode, Eq. (|3.34|) 
remains valid for D > 1 with z^ = z^ = 1. In this case 
the scaling dimension of the purely fermionic four-point 
vertex is £)( 4 >°) = 1 — D and the scaling dimension of the 
three-legged vertex with two fermion legs and one boson 
leg is Z)( 2,1 ) = (1 — D)/2. Both vertices become irrel- 
evant in D > 1. As discussed in the following section, 
this means that the random-phase approximation (RPA) 
for the effective interaction, as well as the so-called GW 
approximation — for the fermionic self-energy, are quali- 
tatively correct in D > 1. 

E. A simple truncation scheme: Keeping only the 
skeleton elements for two-point functions 

In order to solve the flow equations explicitly, one 
is forced to truncate the infinite hierarchy of flow 
equations. In the one-particle irreducible version of 
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(a) 



(b) 





(c) 




T 




FIG. 8: Skeleton diagrams for (a) the one-particle irreducible 
ferrnionic self-energy; (b) the one-interaction-line irreducible 
polarization; and (c) the three-legged vertex with two fermion 
legs and one boson leg. The small black circle denotes the 
bare three-legged vertex. Thin lines denote external legs. The 
other graphical elements are the same as in Fig. 2] 



the purely ferrnionic functional RG it is common 
practicei£ii£iiLi2i2§i2£ to retain only vertices up to the 
four-point vertex and set all higher order vertices equal 
to zero. Our approach offers new possibilities for trun- 
cation schemes. Consider the skeleton graphs^ for 
the one-particle irreducible ferrnionic self-energy and 
the one-interaction-line irreducible polarization shown in 
Figs. Ufa) and (b). The skeleton graphs contain three 
basic elements: the exact ferrnionic Green's function, 
the exact bosonic Green's function (i.e., the effective 
screened interaction), and the three-legged vertex with 
two fermion legs and one boson leg. A systematic deriva- 
tion of the skeleton expansion for the vertices in our cou- 
pled Fermi-Bose theory is presented in Appendix B. One 
advantage of our RG approach (as compared with more 
conventional methods involving only ferrnionic fields) is 
that it yields directly the flow equations for basic ele- 
ments appearing in the skeleton graphs for the self-energy 
and the polarization shown in Fig. Of course, in prin- 
ciple the three-legged vertex can be obtained from the 
vertex with four fermion legs with the help of the skele- 
ton graph shown in Fig. [Sfc). However, calculating the 
three-legged vertex from the four-legged vertex in this 
way involves an intermediate integration, which requires 
the knowledge of the momentum and frequency depen- 
dence of the four-legged vertex. Unfortunately, in prac- 
tice the purely ferrnionic functional RG equations have to 
be severely truncated so that up to now it has not been 
possible to keep track of the frequency dependence of the 
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(b) 
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FIG. 9: Truncation of the flow equations for (a) ferrnionic 
self-energy, (b) irreducible polarization, and (c) three-legged 
vertex which sets all other vertices equal to zero. The internal 
lines are full propagators, which depend on the self-energies 

r< 2 -°) = e and r(°- 2 > = n. 



four-legged fermion vertex within the purely ferrnionic 
functional RG. 

To obtain a closed system of RG equations involving 
only the skeleton elements, let us retain only the vertices 
T, a (K), U a (K) ajidrW(K+Ka;Ka;Ka) on the right- 
hand sides of the exact flow equations for these quantities 
shown in Figs. 00 and0 and set all other vertices equal 
to zero. The resulting closed system of flow equations is 
shown graphically in Fig.0 Explicitly, the flow equations 



F crcr (K)G cr (K + K) + F (J< j{K)G <J {K + K) 

K L 

;r (2,i) ^ K + K(J . Ka)T {2 - 1] (Ka; K + Ka; -Ka) , 

(3.35) 
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a A n CT (if) 



-c 



K 
-(2,1) 



Ga(K)G a (K + K) + Gcr(K)G a (K + K) 

(2,1)/ 



xT^ L '(K + Ka; Ka; Ka)T (z > L> (Ka; K + Ka; -Ka) 



(3.36) 



9 A r (24) (K + Ka; Ka; Ka) 



[ \P aa (K')G a {K + K')G a {K + K + K') 

J K' L 

+F aa (K')G a (K + K')G a (K + K + K') 

+F aa (K')G a (K + K')G a (K + K + K') 

xT^(K + Ka;K + K + K'a; -K'a) 
xT {2 > 1] {K + K + K'a; K + K'a; Ka) 
X T^(K + K'a; Ka; K'a) . (3.37) 

These equations form a closed system of integrodiffer- 
ential equations that can in principle be solved numeri- 
cally. If the initial momentum transfer cutoff Ao is cho- 
sen larger than the maximal momentum transferred by 
the bare interaction, and if the initial bandwidth cutoff 
voAq is larger than the bandwidth of the bare energy 
dispersion, then the initial conditions are Y> a (K) Aa = 0, 
n^l^Ao = 0, and T^iK + Ka;Ka;Ka) Ao = i. A 
numerical solution of these coupled equations seems to 
be a difficult task, which we shall not attempt in this 
work. Note, however, that in Sec. IllllTl we have argued 
that for regular interactions in dimensions D > 1 the 
three-legged vertex is actually irrelevant in the RG sense. 
Hence, we expect that the qualitatively correct behav- 
ior of the fermionic self-energy and of the polarization 
can be obtained by ignoring the flow of the three-legged 
vertex, setting r( 2:1 ) — > i. If we further ignore interac- 
tion corrections to the internal propagators in the flow 
equation (|3.3(il) for the polarization, it is easy to see that 
the solution of this equation is nothing but the nonin- 
teracting polarization. This is equivalent with the RPA 
for the effective interaction. Substituting this into the 
flow equation l|3.35|l for the self-energy and ignoring again 
self-energy corrections to the internal Green's functions, 
we obtain the non-self-consistent GW approximation 56 
for the fermionic self-energy. For regular interactions in 
D > 1 wc therefore expect that the RPA and the GW 
approximation are qualitatively correct. However, for 
strong bare interactions quantitatively accurate results 
can only be expected if the vertex corrections described 
by Eq. I|3.37|l are at least approximately taken into ac- 
count. 

We shall consider this problem again in Sec. IIVDI 
where we discuss truncations of an expansion based on 
relevance. To lowest order, this approximation will agree 
with Eqs. (|3.35l) - (|3.37|l when the dependence of the ver- 
tex r^ 2,1 ) on momenta and frequencies is ignored. There, 
we use the resulting equations to calculate an approx- 
imation to the electronic Green's function of the one- 
dimensional Tomonaga-Luttinger model. Amazingly, 



this simple truncation is sufficient to reproduce the cor- 
rect anomalous dimension known from bosonization even 
for large values of the bare coupling. 



IV. THE MOMENTUM-TRANSFER CUTOFF 
AS FLOW PARAMETER 



The truncation discussed in Sec. IIII El violates the 
Ward identities relating vertices with different numbers 
of external legs (for a self-contained derivation of the 
Ward identities within the framework of our functional 
integral approach; see Appendix Q. Moreover, even if 
we do not truncate the exact hierarchy of flow equations 
shown in Figs. El El and [7| the Ward identities are vi- 
olated for any finite value of the bandwidth-cutoff vqA, 
because the cutoff leads to a violation of the underly- 
ing gauge symmetry. We can thus only expect the Ward 
identities to be restored in the limit vqA — » 0. Recall 
that in the Tomonaga-Luttinger model the Ward identi- 
ties are valid in the strict sense only in the presence of 
the Dirac sea, implying that the ultraviolet cutoff t>oAo 
has been removed. The Ward identities and the underly- 
ing asymptotic conservation laws are crucial for the ex- 
act solubility of the TLM 43,44 and its higher-dimensional 
generalizationi 4 ^ 4 ^ 5 ^^ In order to reproduce the exact 
solution of the TLM known from bosonization within the 
functional RG, it is very important to have RG flow equa- 
tions which are consistent with the Ward identities, even 
for finite values of the cutoff. In this section we show 
that this requirement is fulfilled if we work in our mixed 
Fermi-Bose RG with a momentum-transfer cutoff A only 
and take the limit vqA — * of a vanishing bandwidth 
cutoff. 



Exact flow equations for momentum-transfer 
cutoff 



As already briefly mentioned at the end of Sec. IIII Gl if 
we work with a momentum transfer cutoff A only, then all 
diagrams with a slash on an internal fermionic Green's 
function [corresponding to the fermionic component of 
the single-scale propagator given in Eq. I|3.13|l ] on the 
right-hand sides of the exact flow equations shown in 
Figs. El El and[71should be omitted. The exact flow equa- 
tions for the electronic self-energy and the irreducible po- 
larization then reduce to 

d A X a (K) = i / F aa (K)rW{Ka;Ka;Ka,-Ka) 

F aa (K)G a {K + K)T {2 ^ (K + Ka; Ka; Ka) 

xT^(Ka; K + Ka; -Ka) , (4.1) 





FIG. 10: Exact flow equations for (a) the fermionic self- 
energy and (b) the irreducible polarization in the momentum- 
transfer cutoff scheme. 





FIG. 11: Exact flow equations for the three-legged vertex 
with two fermion legs and one boson leg in the momentum- 
transfer cutoff scheme. 



z JK' 

- [ F aa (K')F aa (K + K')T^ 3 \-Ka,K + K'a,-K'a) 

JK' 

xT i0 ^(K'a, -K - K'a, Ka) . (4.2) 
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FIG. 12: Initial condition for the pure boson vertices in the 
momentum transfer cutoff scheme. The sum is taken over 
the m\ permutations of the labels of the external legs. For 
linearized energy dispersion, all symmetrized closed fermion 
loops with more than two external legs vanish. 

These equations are shown graphically in Fig. 111)1 A 
graphical representation of the corresponding exact flow 
equation for the three-legged vertex is shown in Fig. 1111 
Still, these flow equations look rather complicated. Since 
we have imposed a cutoff only in the momentum trans- 
fered by the bosons, the initial conditions at scale Ao are 
nontrivial. The initial value of the three-legged vertex is 
still r^' 1 ' = i, but the pure boson vertices r' '" 1 ) with m 
external legs are initially given by the symmetrized closed 
fermion loops shown in Fig.^| All other vertices vanish 
at the initial scale Ao. An essential simplification occurs 
now if we linearize the energy dispersion relative to the 
Fermi surface. If the initial momentum transfer cutoff 
Ao is small compared with the typical Fermi momentum, 
then we may set all pure boson vertices r' '" 1 ' with more 
than two external boson legs (m > 2) equal to zero. This 
is nothing but the closed loop theorem^&iM&iiliffi which 
is valid exactly for the one-dimensional TLM (where the 
energy dispersion is linear by definition). In higher di- 
mensions, the closed loop theorem is valid to a very good 
approximation as long as the linearization of the energy 
dispersion is justified within a given sectorization of the 
Fermi surface and scattering processes that transfer mo- 
mentum between different sectors of the Fermi surface 
can be neglectedi 47 ' 49 Note that the closed loop theorem 
is consistent with the momentum-transfer cutoff flow, be- 
cause pure boson vertices r^ '™- 1 with m > 3 are not 
generated if they initially vanish. 

Assuming the validity of the closed loop theorem, the 
right-hand side of the flow equation l|4.2fl for the polar- 
ization vanishes identically, because it depends only on 
boson vertices with more than two external legs. Phys- 
ically, this means that there are no corrections to the 
nonintcracting polarization, so that the RPA for the ef- 
fective interaction is exact. This is of course well known 
since the pioneering work by Dzyaloshinskii and Larkin. 43 
Moreover, the last three diagrams in the flow equation for 
pC 2 ' 1 ) shown in Fig. ^] also vanish, because they contain 
the vertex r*- 0,3 ' . However, the remaining diagrams in 
Fig. HUT a - ) and Fig. ^2 still look quite complicated, so 
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that we still have to solve an infinite hierarchy of cou- 
pled flow equations. In the next subsection we show how 
this infinite system of coupled integrodifferential equa- 
tions can be solved exactly. 

B. Ward identities as solutions of the infinite 
hierarchy of flow equations 

Let us consider the terms on the right-hand sides of the 
flow equations for the vertices I^ 2 '™) with two external 



fermion legs and an arbitrary number of boson legs. As- 
suming the validity of the closed loop theorem, all pure 
boson vertices r(° ,m ) with m > 2 vanish. From Fig. llOf a) 
and Fig. ^2 it is clear that in general the right-hand side 
of the flow equation for ^Ar^ 2 '™^ depends on r( 2,m+2 ) 
and on all T^ 2,m ) with m' < m. In fact, from our general 
expression for the flow of the totally symmetrized vertices 
given in Eq. 13.16[1 . we can derive the flow equations for 
the vertices r^ 2 '" 1 ) with arbitrary m in closed form (we 
omit for simplicity the degeneracy index a), 



1 f °° 00 d 

dA r( 2 <™)(K>;K;K u ...,K m ) = ~ F aa {K)T^ m+2 \K' ; K; -K, K, K u . . . , K m ) + V V 

JK 1=2 mi mi=l lu 

f F(K)T^ m ^ (k 1 ; Ki;K p{1) , . . . , K P(mi) , -I<) G^)^ 2 ^ fa K 2 ; K P{mi+1) , . . . , K P{mi+m2 

p J K 

xG(K 2 ) • . . . • G(Jr,_i)r( 2 ' m ' +1 > (Ki-i;K; K, K P{m _ mi+1) , . . . , K P(rn) ) , (4.3) 



r 



where we have defined 



K' = K + J2 K i 



mi+...+mi 

Ki = K' + K- R p(j) - 

i=l j=l 

(4.4) 

and P denotes a permutation of {1, ... , m}. A graphical 
representation of Eq. 14.3fl is shown in Fig. ^| Note that 
the flow equation (|4.1|) for the irreducible self-energy is a 
special case of Eq. (|4.3f) for m = 0. 

We are now facing the problem of solving the infinite 
hierarchy of coupled flow equations given by Eq. 14.311 . 
In view of the fact that these equations are exact and 
that in one dimension the single-particle Green's function 
of the TLM can be calculated exactly via bosonization, 
we expect that this infinite hierarchy of flow equations 
can also be solved exactly. Indeed, the solutions of these 
equations are nothing but infinitely many Ward identities 
relating the vertex r( 2,m ) with two fermion legs and m 
boson legs to the vertex r^ 2 '™" 1 ) with one boson leg less. 
Wc derive these Ward identities within the framework 



of our functional integral approach in Appendix C. For 
rn = 1 the Ward identity is well knowni244iiS4i 

G(K + K)T^ 1 \K + K;K;K)G(K) = 



G(K + K)-G(K) . (4.5) 



ilu — v Pa ■ k 



Here v P<a is the Fermi velocity associated with the inde- 
pendent fermionic label K = (k, iuj), where |k — ki? CT | <C 
IkF.oJ. The Ward identity (|4.5|l has been used in Refs.^3 
and |45| to close the skeleton equation for the self-energy 
and thus obtain the exact Green's function of the TLM 
without invoking the machinery of bosonization. A Ward 
identity for r^ 4,1 ) has also been used to prove the van- 
ishing of the renormalization group f3 function for the 
TLM. 46 However, for solving the TLM exactly within 
the framework of the functional RG, we need the Ward 
identities for all vertices r^ 2 '™) with m > 1. As shown in 
Appendix C, for linear energy dispersion we have 



T^(K';K-K 1 ,... 1 K rn ) = 



iuji - vf,o- ■ k/ 

_p(2,m 



r (2,m-i) ^ Rl . K + k v ,K u ..., K^ u K l+1 , ...,K, 
Q (K 1 - K i; K; K u ..., Ki-!,Ki +1 , ...,K 



(4.6) 
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2,m 




FIG. 13: Diagrammatic representation of the flow equation l|4.^ of vertices with two fermion legs and a general number of 
boson legs provided the pure boson vertices with more than two external legs vanish, as implied by the closed loop theorem. 



where 1 < I < m. For clarity let us write down here the special case m = 2, 



r (2,2) f K + Kl + x 2 . K . K K \ = 



iu)\ - vp a • ki 



-(2,1) 



K + Ki + K 2 ; K + Kr,K 2 ) - T^ 2 ^ [K + K 2 ; K; K 2 



iu 2 - VP a ■ k 2 



r^ 2 ' 1 ) (K + K 1 + K 2 :K + K 2 - - (k + K i; K; K x 



(4.7) 



Diagrammatic representations of the Ward identities 
given in Eqs. I|4.5|l and l|4.6|l are shown in Fig. To 
prove that these Ward identities indeed solve our infinite 
system of flow equations given by Eqs. gdj) and 
we start from the flow equation for r( 2 < m+1 ). Substitut- 
ing on both sides of this exact flow equation the Ward 
identities, we can reduce it to a new flow equation in- 
volving only vertices where the number of boson legs is 
reduced by one, but with an external bosonic momentum 
entering the vertices at various places. Graphically, we 
indicate the place where the bosonic momentum enters 
the vertex by a double slash, as shown in Fig. ^] The 
important point is now that all diagrams with double 
slashes attached to intermediate Green's functions cancel 
due to the fact that all vertices r( 2,m ) can be expressed 
in terms of a difference of vertices r^ 2 '" 1-1 ), with a same 
prefactor that is independent of m. Graphically, only the 




FIG. 14: (a) Diagrammatic representation of the Ward iden- 
tity 14.51 for the three-legged vertex and (b) of the Ward 
identity 14.611 for the vertex with two fermion legs and m > 1 
boson legs. The small arrow indicates the place in the dia- 
gram where the external bosonic energy-momentum enters. A 
double slash to the right of an arrow means that the bosonic 
momentum is added before the corresponding Green's func- 
tion, while a double slash to the left of an arrow means that 
the momentum is added after the Green's function. 



diagrams with a double slash attached to the leftmost or 
rightmost Green's function survive. Canceling the com- 
mon prefactor, it is then easy to see that the RG equation 
derived in this way from the functional RG equation for 
r (2,m+i) ig noting but the exact RG equation for r( 2 < m ). 
Hence, the Ward identities provide relations between the 
vertices r( 2,m ) that are consistent with the relations im- 
plied by the exact hierarchy of RG flow equations in the 
momentum-transfer cutoff scheme. In other words, the 
Ward identities are the solutions of the infinite hierarchy 
of flow equations! 

C. Exact solution of the Tomonaga-Luttinger 
model via the exact RG 

Given the cascade of Ward identities 14.5fl and l|4.6|) 
we can close the integrodifferential equation 14.1fl for the 
irreducible self-energy. Note that this equation involves 
both the three-legged vertex and the four-legged vertex 
with two fermion legs and two boson legs, so that the 
Ward identity l|4.5[) is not sufficient to close the flow equa- 
tion. Of course, if one is only interested in calculating 
the Green's function of the TLM, it is simpler to start 
from the skeleton equation for the self-energy shown in 
Fig. |SJ which can be closed by means of the Ward iden- 
tity 1)4. 5fl for the three-legged vertex only. Nevertheless, 
it is instructive to see how the exact solution emerges 
within the framework of the functional RG. Substituting 
Eqs. 14.5f> and l|4.7[) into Eq. I|4.1|l . we obtain the following 
integrodifferential equation for the electronic self-energy: 

d^ a {K) = G?{K) I r - FaAK \, 2 
Jk (iu - Vf,<t ■ k) 2 

x [G„{K) - G„(K + Kj\ . (4.8) 

Here the index a labels not only the different spin 
species, but also the different patches of the sector- 
ized Fermi surfaced For example, for the spinless case 
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a = ±kp. Using the fact that in the momentum-transfer 
cutoff scheme G 2 9aS = d\G we can alternatively write 
Eq. H4.8J1 as a linear integrodifferential equation for the 
fermionic Green's function, 



d A G a (K) = f 

Jb 



K (W - VF,o- ■ k) 2 

x [G a {K) - G a (K + K)} . (4.9) 

If we had simply set the vertex I^ 2,2 ) equal to zero in 
Eq (|4.1|l and had then closed this equation by means of 
the Ward identity l]4.5[l . we would have obtained a non- 
linear equation. Thus, the linearity of Eq. I|4.9J) is the 
result of a cancellation of nonlinear terms arising from 
both Ward identities (|4.5|l and (|4.7|) . Because the second 
term on the right-hand side of Eq. 14.9fl is a convolution, 
we can easily solve this equation by means of a Fourier 
transformation to imaginary time and real space. Defin- 
ing 



G a {X) = [ e l ^ r -^G a (K) 
Jk 



(4.10) 



e-&- r - OT K. F ° ai ~ K) , (4.11) 



A' 



lUJ - Vjr i<7 ■ k) 2 



where X — (r, r), the flow equation (|4.9|l is transformed 
to 



[d A + H A , a (X) - H\ ;f j (0)J G a {X)=0. 
This implies the conservation law 

J\a' [H K ,. a {X) - ffAv(O)] j G a (X) 





exp < 


I 









(4.12) 

= 0. 

(4.13) 



Integrating from A = to A = Aq, we obtain 

G a {X) = GoAX) exp [Qa(X)} , (4.14) 



with 



and 



Q (T (X) = S !T (0)-S a (X), 



(4.15) 



S a (X) 



dMHv, a {X) 
9(Ao - \k\)F aa (K) 



K («W - VF,cr • k) 2 



cos(k ■ r — lot) , (4.16) 



where we have used the invariance of the RPA interac- 
tion F(K) under K — > —K. The solution in Eqs. I|4.14fl - 
(|4.1(il) is well known from the functional integral approach 
to bosonizationiS^SiSiffi where Q a (X) arises as a Debye- 
Waller factor from Gaussian averaging over the distribu- 
tion of the Hubbard-Stratonovich field. In one dimension, 
Eqs. (|4.14|) - (|4.16|) can be shownii to be equivalent to the 



exact solution for the Green's function of the Tomonaga- 
Luttingcr model obtained via conventional bosonization. 

Once the exact single-particle Green's function is 
known, the Ward identities in Eqs. I|4.5|) and (|4.6|l itera- 
tively yield expressions for the vertices I^ 2 '™' that solve 
the whole hierarchy of flow equations for the vertices with 
two fermion and an arbitrary number of boson legs. In 
principle, the method described in this section can be 
applied also to vertices with more than two fermion legs. 
For example, the right-hand sides of the flow equations 
for the vertices r( 4,m ) contain only vertices with no more 
than four fermion legs. Ward identities for these vertices 
would again yield a solution of this complete hierarchy, 
once the vertices r' 2 '™) are known. This procedure can 
be iterated to obtain vertices with an arbitrary number 
of external legs using at each step the complete flow of 
vertices with two fewer fermion legs obtained in the pre- 
vious step. We have thus devised a method to obtain 
all correlation functions of the TLM entirely within the 
framework of the functional RG. 



D. Truncation scheme based on relevance 

The structure of the exact Green's function of the 
TLM and the corresponding spectral function A(k,ui) = 
—ir~ 1 ImG(k,uj + iO) depend crucially on the Ward iden- 
tities discussed above, which in turn are only valid if the 
energy dispersion is strictly linear. In order to assess the 
validity of the linearization of the energy dispersion, it 
is important to develop truncations of the exact hierar- 
chy of flow equations that do not explicitly make use of 
the validity of the asymptotic Ward identities. We now 
propose such a truncation scheme. 

The coefficients generated in the expansion of a given 
vertex T^ m \K' 1: . . . , K' n ; K u . . . , K n ; K u . . . , K m ) in 
powers of frequencies and momenta have decreasing scal- 
ing dimensions, so that the most relevant part of any ver- 
tex is obtained by setting all momenta and frequencies 
equal to zero. This classification leads to a simple trun- 
cation scheme: We retain only those vertices whose lead- 
ing (momentum- and frequency-independent) part has a 
positive or vanishing scaling dimension, corresponding to 
relevant or marginal couplings in the usual RG jargon. 
In the context of calculating the critical temperature of 
the weakly interacting Bose gas in three dimensions, such 
a truncation procedure has recently been shown to give 
very accurate results^ 

To begin with, let us classify all couplings according to 
their relevance. With the rescaling defined in Sec. IIII Dl 
for D = = z v = 1, the scaling dimensions of the 
vertices f( 2n - m ) are £>( 2 ™. m ) = 2 - n - m, see Eq. f3~31) . 
Hence the vertex r < - 2,2 - 1 as well as the vertices f^ ' 3 * 1 and 
p(o>4) 5 wnose unrescaled versions appear on the right- 
hand sides of Eqs. I|4.1|l and ^4.2(1 . are irrelevant in the 
RG sense. In contrast, the momentum- and frequency- 
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independent part of the three-legged vertex, 
/ A 



7, =f( 2 - 1 )(0;0;0) 



Z l Y^ 1 \K F ;K F ;Q) , 



(4.17) 

is marginal^ Here Kp = (±kp,w = 0). From the gen- 
eral flow equations l|3.26[l and (|3.30|) for the rescaled ver- 
tices we obtain the following exact flow equation for the 
rescaled self-energy defined in Eq. i|3.24[l : 



r(2,o) 



(Q) 



with [see Eq. pi]) ] 



r 



(2:0) 



(Q) 



Ad A T^ 0) (K) 



(4.18) 



(4.19) 



We restrict ourselves to spinless fermions here and choose 
fl\ = — wpA, so that with vq = (ttvf)~ 1 the prefac- 
tor in Eq. (|4"TT|) turns out to be (^^) 1/2 = tt 1/2 . As 
usual, the fermionic wave- function renormalization factor 
Zi is defined via 



1 



d(iu>) 



K=0 



= 1 



8(ie) 



(4.20) 



Q=0 



According to Eq. I|3.25|) the wave-function renormaliza- 
tion Zi satisfies the flow equation 



diZi = -rjiZi , 



(4.21) 



where the flowing anomalous dimension of the fermion 
fields is given by 



dt\ 2fi) (Q) 



d(ie) 



(4.22) 



Q=0 



According to Eq. (|4. the constant part of the self- 
energy, 

fi = ti(0), (4.23) 

is relevant and satisfies 

9m = (l-J70n+ff' O) (O). (4.24) 

In general, f/ will only flow into the fixed point if the 
initial coupling f o is properly fine-tuned. Apart from Zi , 
there are two more marginal couplings. The first is the 
Fermi velocity renormalization factor— 



z t + 



dt^Q) 



dq 



(4.25) 



3=0 



and the second marginal coupling is the momentum- 
and frequency-independent part ji of the rescaled three- 
legged vertex given in Eq. I|4.17|l . The exact flow equa- 
tions for vi and 7; are 



diiii = —TjlVl + 



afp' 0) (Q) 



dq 



(4.26) 



Q=0 



and 



-(2,1: 



(0;0;0) . 



(4.27) 



If we retain only relevant and marginal couplings, then 
in the momentum-transfer cutoff scheme the rescaled 
fermionic Green's function defined in Eq. (|3.32|1 is in 
D = 1 simply approximated by 



G(Q) 



1 



ie - vtq - n 



(4.28) 



In order to make progress, we have to approximate the 

inhomogeneities ff ' 0) (Q) and f | 2,1) (0; 0; 0). In SecHHH 
we have proposed an approximation scheme which retains 
only the skeleton elements of the two-point functions. In 
the momentum-transfer cutoff scheme, the corresponding 

1 13.37(1 further simplify because 



flow equations 
we should omit all terms involving the fermionic single- 
scale propagator. Unfortunately, the resulting nonlinear 
integrodifferential equations still cannot be solved ana- 
lytically. In order to simplify these equations further, let 
us replace the three-legged vertex on the right-hand sides 
of these equations by its marginal part. In this approxi- 
mation we obtain, from Eq. (|3.35() . 



f{ 2,0) (Q)»7? / F(Q)G(Q + Q) , 
Jq 

and from Eq. Q3.37[) . 

r{ a ^(0;0;0)w7f / P(Q)G 2 (Q) . 



(4.29) 



(4.30) 



In order to be consistent, we should approximate G(Q) 
in Eqs. flO^ and by Eq. g^HJ. Then it is easy 

to see that the second term on the right-hand sides of 
the flow equations l|4.26|l and (|4.27|) exactly cancels the 
contribution from the anomalous dimension, so that 



dai = , dm = 



(4.31) 



For explicit calculations, let us assume that the usual 
couplings of the TLM 52 are 52 = ff4 = /o, so that 



F(Q) = 6(1 - \q\) 



fo(f 



(1 + f Q )q 2 + e 2 



(4.32) 



where f = Vofo- From Eqs. (j4.22jl and (j4.29jl we then 
find that the anomalous dimension 77 = 77/ does not flow 
and is given by££ 



/o 2 



2\/l + /o 



1 + fo + 1 



(4.33) 



which agrees exactly with the bosonization result 4L£i 
We emphasize that Eq. I|4.33[) is the correct anomalous 
dimension of the TLM, even for Jo 3> 1, so that, at least 
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as far as the calculation of rj is concerned, the validity of 
our simple truncation is not restricted to the weak cou- 
pling regime. Recall that the restriction to weak coupling 
is one of the shortcomings of the conventional fermionic 
functional R G , 15 i 16 , 17 i 18 , 20 i 21 , 22 i 23 , 24 i 25 i 26 i 27 which was 
implemented for the TLM in Ref. 0. Because rj 
is finite, the running vertex I^ 2,1 ) (Kp ; Kp; 0) without 
wave-function renormalization actually diverges for A — > 
0. However, the properly renormalized vertex 7; oc 
ZiT 1 - 2 -' 1 ^ (Kp; Kp; 0) remains finite due to the vanishing 
wave-function renormalization 

^ = e-*=(£f (4-34) 

for I — > 00. Integrating the flow equation l|4.18|l for 
the self-energy with the inhomogeneity approximated by 
Eqs. 14.29JI and (|4.28(l . we obtain, after going back to 
physical variables 63 

£(fcF + k, iuo) — 

r Aa dfc c°° du)_ /A y / RPA (fc,^) 

J_ Ao 2ir J_ 00 2ir \\k\J i(u + Q) - Vp(k + k) ' 

(4.35) 

where the RPA screened interaction is 

r\kM = ■ (4.36) 



Here v c = vpyl + fo is the velocity of collective 
charge excitations. Equation l|4.35[) resembles the 
GW approximation^ but with the RPA interaction 
multiplied by an additional singular vertex correction 
(Ao/\k\)' n . The explicit evaluation of Eq. (|4.35|) is rather 
tedious and will not be further discussed in this work. 
The resulting spectral function A(k,cu) agrees at k = kp 
with the bosonization result (even at strong coupling), 
but has the wrong threshold singularities for \to\ — > 
v c \(k ± kp)\. So far we have not been able to find a 
reasonably simple truncation of the exact flow equations 
which completely produces the spectral line shape of 
A(k,ui), as predicted by bosonization or by our exact so- 
lution presented in the previous section. Whether a self- 
consistent numerical solution of the truncation discussed 
in Sec. lIII El [see Eqs. (|3~35|) - (|3~i?7|) ] would reproduce the 
correct spectral line shape or not remains an open prob- 
lem. The numerical solution of these equations seems to 
be rather difficult and is beyond the scope of this work. 

V. SUMMARY AND OUTLOOK 

In this work we have developed a new formulation of 
the functional RG for interacting fcrmions, which is based 
on the explicit introduction of collective bosonic degrees 
of freedom via a suitable Hubbard-Stratonovich transfor- 
mation. A similar strategy has been used previously in 



Refs.l3^.l38l and 13^. However, on the technical level the 
practical implementation of this method presented here 
differs considerable from previous works. We have payed 
particular attention to asymptotic Ward identities, which 
play a crucial role if the interaction is dominated by small 
momentum transfers. In one dimension, this is the key 
to obtain the exact solution of the Tomonaga-Luttinger 
model entirely within the functional RG. By using the 
momentum transfer associated with the bosonic field as 
a cutoff parameter, we have formulated the functional 
RG in such a way that the RG flow does not violate the 
Ward identities. In fact, we have shown that Ward iden- 
tities emerge as the solution of the infinite hierarchy of 
coupled RG flow equations for the one-line irreducible 
vertices involving two external fermion legs and an arbi- 
trary number of boson legs. In principle this method can 
be iterated to obtain all correlation functions of the TLM 
entirely within the framework of the functional RG. 

Here we have mainly laid the theoretical foundation of 
our approach and developed an efficient method to keep 
track of all terms. In future work, we are planning to 
apply our technique to other physically interesting prob- 
lems. Let us mention some problems where it might be 
advantageous to use our approach: 

(a) Strong coupling fixed points. One of the big draw- 
backs of the conventional (purely fermionic) functional 
RG used by many authora 15 i 16 i 17 i 18 i 23 i 24 i 25 i 26 i 27 is that in 
practice the frequency dependence of the four-point ver- 
tex r^ 4 ) has to be neglected, so that the wave-function 
renormalization factor is Z = 1. Although the resulting 
runaway flow of the vertices to strong coupling at a finite 
scale can be interpreted in terms of corresponding insta- 
bilities, there is the possibility that for small Z the renor- 
malized effective interaction Z 2 T^ remains finite even 
though the vertex r^ 4 -* without wave-function renormal- 
ization seems to diverged In our approach, the effective 
interaction acquires a frequency dependence, even within 
the lowest-order approximation. In fact, if we ignore ver- 
tex corrections, the effective interaction is simply given 
by the RPA. Hence, strong coupling fixed points might 
be accessible within our approach. Recall that the rather 
simple truncation of Sec. HVDl gave the exact anomalous 
dimension of the TLM for arbitrary strength of the in- 
teraction. Possibly, more elaborate truncations of the 
exact hierarchy of RG flow equations (for example, the 
truncation based on retaining skeleton elements of the 
two-point functions discussed in Sec. IIII El see Fig. |3J) 
will give accurate results for the spectral properties. 

(b) Nonuniversal effects in one- dimensional metals. If 
we do not linearize the energy dispersion in one dimen- 
sion, there should be a finite momentum scale k c (de- 
pending on the interaction and the band curvature) be- 
low which typical scaling behavior predicted by the TLM 
emerges. The calculation of k c as well as the associ- 
ated nonuniversal spectral line shape are difficult within 
bosonization^ On the other hand, within the framework 
of the functional RG the inclusion of irrelevant couplings 
is certainly possible, so that with our method it might be 
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possible to shed some new light onto this old problem. 
For an explicit calculation of an entire crossover scal- 
ing function between the critical regime and the short- 
wavelength regime of interacting bosons in D = 3, see 
Ref. |^. An analogous calculation of the dynamic scal- 
ing functions for interacting fermions in one dimension 
remains to be done. 

(c) Itinerant ferromagnetism. Spontaneous ferromag- 
netism in Fermi systems is driven by sufficiently strong 
interactions involving small momentum transfers. As- 
suming a given form of the ferromagnetic susceptibility, 
Altshuler, Ioffe, and Millis 66 concluded on the basis of 
an elaborate diagrammatic analysis that in the vicin- 
ity of the paramagnetic-ferromagnetic quantum-critical 
point in dimensions D — 2 a simple one-loop calcula- 
tion already yields the correct qualitative behavior of the 
electronic self-energy. If this is true, then in this problem 
vertex corrections are irrelevant. Unfortunately, due to 
the peculiar momentum and frequency dependence of the 
ferromagnetic susceptibility x(k,a)) at the quantum crit- 
ical point, the assumption of asymptotic velocity conser- 
vation [see Eq. I|C13|) in Appendix C] leading to the Ward 
identity 14.5ll - H4.7l) is not justified. Note, however, that 
in Ref. |66| the form of the susceptibility is assumed to be 
given. The feedback of the collective ferromagnetic fluc- 
tuations on the non-Fermi-liquid form of the electronic 
properties has not been discussed. The fact that in D < 3 
the leading interaction corrections to the inverse sus- 
ceptibility x -1 (k,d>) generate a nonanalytic momentum 
dependenc e 67 ! 68 suggests that the problem should be re- 
considered taking the interplay between fermionic single- 
particle excitations and collective magnetic fluctuations 
into account. The formalism developed in this work 
might be suitable to shed some new light also onto this 
problem. The electronic properties of a two-dimensional 
Fermi system in the vicinity of a ferromagnetic instability 
have recently been studied in Refill However, these au- 
thors focused on the finite-temperature properties of the 
phase transition; they did not attempt to calculate the 
fermionic single-particle Green's function in the vicinity 
of the zero-temperature phase transition. Note also that 
for sufficiently strong interactions even one-dimensional 
fermions can in principle have a ferromagnetic instability 
if the energy dispersion is nonlinear 

(d) Quantum phase transitions and symmetry breaking. 
Our method unifies the traditional approach to quantum- 
phase transitions pioneered by HertaS with the modern 
developments in the field of fermionic functional RG, so 
that it might simplify the theoretical description of quan- 
tum phase transitions in situations where the fermions 
cannot be completely integrated out. In order to describe 
quantum phase transitions in interacting Fermi sys- 
tems within the framework of the traditional Ginzburg- 
Landau- Wilson approach, all soft modes in the system 
should be explicitly retained^ In the purely fermionic 
functional RGiVVVVVO, 21 ! 22 ,^ 24 ! 2 ^ 6 , 2 ? symme try 
breaking manifests itself via the divergence of the rel- 
evant order-parameter susceptibility; the symmetry bro- 



ken phase is difficult to describe within this approach. 
On the other hand, in our approach the order parame- 
ter can be introduced explicitly as a bosonic field, which 
acquires a vacuum expectation value in the symmetry 
broken phase. Previously, a similar approach has been 
developed in Ref. |3^ to study antiferromagnetism in the 
two-dimensional Hubbard model. 

In summary, we believe that the formulation of the 
exact functional RG presented in this work will be quite 
useful in many different physical contexts. 
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APPENDIX A: TREE EXPANSION OF 
CONNECTED GREEN'S FUNCTIONS IN TERMS 
OF ONE-LINE IRREDUCIBLE VERTICES 

In this Appendix we show explicitly that the vertices 
r^...Q. n defined in terms of the functional Taylor expan- 
sion of r[$] in Eq. 12. 4411 are indeed one-line irreducible. 
This is usually 10 done graphically by taking higher-order 
derivatives of the relation (|2.4(J|) between the second func- 
tional derivatives of £[<&] and Q C {J]. With the help of our 
compact notation we can even give the tree expansion of 
the connected Green's function in terms of one-line irre- 
ducible vertices in closed form. To do so, it is advanta- 
geous to define the functional 



U = 









T 




T 




5<P5<P 


■5=0- 









(Al) 



which is a matrix in superindex space. With this defini- 
tion, we have 



5<f>5<S> 



V 1 



[G-r , 



(A2) 



so that 



-G T [1-V T G T }~ 1 

OO 

-^TG T (U T G T ) ; . (A3) 



1=0 
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From Eq. I|2.40[l we then obtain 



SJ6J 



s&>c 



ZG T [1-U T G T ]- 1 



5$5<S> 
]TZG T (U T G 



(A4) 



We now expand both sides of Eq. (| A4I> in powers of the 
sources J and compare coefficients. For the matrix on 
the left-hand side we obtain, from Eq. I|2.22|l . 



5^ 7T\ 



SJSJ ^ m Jtii 

n— 1 



p(n+2) 



where the matrix Gi^.Q^ is defined by 

rp(n+2) i / _ n{n+2) 



(A5) 
(A6) 



On the right-hand side we use Eqs. (|A1(I and (|2.44|) to 
write 

U = E^ / •••/ r(«+ 2 ) an $ ai .....$ an , (A7) 

n=l Jq i Jq " 



On the right-hand side of this rather cumbersome ex- 
pression, only connected correlation functions with a de- 
gree smaller than on the left-hand side appear. One 
can therefore recursively express all connected correla- 
tion functions via their one-line irreducible counterparts. 
Only a finite number of terms contribute on the right- 
hand side. The operator S symmetrizes the expression 
in curly brackets with respect to indices on different cor- 
relation functions, i.e., it generates all permutations of 
the indices with appropriate signs, counting expressions 
only once that are generated by permutations of indices 
on the same vertex. More precisely the action of S is 
given by (to = Ei=i m 



^Qi ,,..,<x m , 1 ;...;a m _ TO ,_|_i,... ) Q! m {-^ai,.. -,at m J 



n 



(Afl) 



where P denotes a permutation of {1, . . . , to} and sgn^ 
is the sign created by permuting field variables according 
to the permutation P, i.e., 



where 



rWn+2) 1 pira-t-a; 



(n+2) 



(A8) 



To compare terms with the same powers of the sources J 
on both sides of Eq. l|A4j) . we need to express the fields 
$ a on the right-hand side of Eq. (|A7|) in terms of the 
sources, using Eqs. I|2.22[) and (|2.35l) . 
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m— 



r (m+l) j T 
y c ,aP 1 ...(3 m J 0i ■■■ J Pr, 



(A9) 

Substituting Eqs. fKEfy , (|A7|) . and (fA"9|) into Eq. (jA"4|) and 

comparing terms with the same powers of the sources 
(after symmetrization), we obtain a general relation be- 
tween the connected and the one-line irreducible correla- 
tion functions, 



,(n+2) 
"c,0i,...,p. 

X 



_ V V - 

t—i t—i nil-.. 

1—0 ni,...,n^ — 1 

Xati. •■■Sal ■ ■ ■ iaL 



E 



m, ,....m' 



ZG T r (ril+2)T g t • • G T r ( " i+2)T g t 



a ,a' 



XiS, 



01 



(m l +1) 



.}■ 



sgn c (P)$ 



(A12) 



A diagrammatic representation of the first few terms of 
the tree expansion generated by Eq. IjAl Ofl is given in 
Fig. 1151 Let us give the corresponding analytic expres- 
sions: If wc set n = in Eq. lAlOfl . then only the term 
with I = contributes, and we obtain 



g?) 



ZG 



G J 



(A13) 

which is Eq. (|2.24|) in matrix form. For n = 1 the single 
term with I = 1, rii = 1, to} — 1 contributes on the 
right-hand side of Eq. (jAlOp . Using ZG = G T the tree 
expansion of the connected Green's function with three 
external legs can be written as 



[G]/3 iai [G] / 3 2Q2 [G]p 3a3 T^ 



c (3) 

(A14) 

Finally, consider the connected Green's function with 
four external legs, corresponding to n = 2 in Eq. ljA10|l . 
In this case the following three terms contribute: 



mi 


term 


/ 


Hi 


TO} 




1.) 


1 


n± = 1 


to} = 2 




2-) 


1 


m = 2 


to} = ra\ = 1 


(A10) 


3.) 


2 


n\ — n 2 = 1 


to} = m\ = 1 
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The corresponding analytic expression is 




FIG. 15: Graphical representation of the relation between 
connected Green's functions and one-line irreducible vertices 
up to the four-point functions. The irreducible vertices are 
represented by shaded oriented circles with the appropriate 
number of legs; see Fig. Q The connected Green's functions 
are drawn as empty oriented circles with a number indicating 
the number of external legs. 
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[G]f3 iai [G}j3 2 a 2 [G]f3 3 a 3 [G]0 Aa4 Taia a aaon 
[G]/3 iai [G]fa aa [GJ/33Q3 [G],g 4a4 

T (3) [Gl r (3) 

aiQ2CK5 L J a 5 a 6 Q6Q3Q4 

S P 3 ;Pi { [G]/3 lQ1 [G] 02 a 2 [G]f3 aaa [G]p i(Xi 
'^aict5Ct4 [G]a 5 a 6 r„g Q2Q3 | . (A15) 



APPENDIX B: DYSON-SCHWINGER 
EQUATIONS AND SKELETON DIAGRAMS 

In this appendix we show how the skeleton diagrams 
for the two-point functions and the three-legged vertex in 
Fig.[S]can be formally derived from the Dyson-Schwinger 
equations of motion. Although the skeleton graphs are 
usually written down directly from topological consid- 
erations of the structure of diagrammatic perturbation 
theory^ it is instructive to see how the skeleton expan- 
sion of the irreducible vertices can be derived formally 
within our functional integral approach. 

The invariance of the generating functional Q[J] of the 
Green's functions defined in Eq. (|2.18l) with respect to 
infinitesimal shifts in the integration variables $ Q implies 
the Dyson-Schwinger equations of motion^ 



Ca J a 



ss 



5Ja 



G[Ja] = . 



(Bl) 



For our coupled Fermi-Bose system with Euclidean action 
S[tp, ip, tp] given by Eqs. (|221 EH EH EE1> involving three 
types of fields, Eq. (jBlfl is actually equivalent with the 
following three equations: 



< 



S&>Q 



6 



= 0, 



(B2) 




= 



(B3) 



k 8] K _ Ka 5J_ Ka 



= 



(B4) 



Expressing these equations in terms of the generating 
functionals G c [j,j,J] of the connected Green's functions 
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and the corresponding generating functional T[ip, ip, <p] of 
the irreducible vertices defined in Eq. 1)2 .42 1) . we obtain 
the Dyson-Schwinger equations of motion in the following 
form: 



ST 
6T 
ST 



K 



K 



K 



4> 



K+K,aYKa 



5]Ka5j K +K,a 

s&>g c 



$JK+K,a&J-K<r 



^K-K.afRa + 



&3K-K,cj&J-Ka 



= o, 

(B5) 
= 0, 
(B6) 
= . 
(B7) 



The second functional derivatives of Q c can be expressed 
in terms of the irreducible vertices using Eq. (53). Tak- 
ing derivatives of Eqs. I|B5HB7|I with respect to the fields 
and then setting the fields equal to zero we obtain the 
desired skeleton expansions of the irreducible vertices. 
Let us start with the skeleton diagram for the self-energy 
shown in Fig.|Hl[a). To derive this, we simply differentiate 
Eq. I|B7|I with respect to ipx'o- Using 



SlpK'aStpKa 



(B8) 



ficlds=0 



we obtain 



5 KK ,V a {K) = i f 
J i\ 



ficlds=0 



(B9) 

From the I = 1 term in the expansion (|A4J) it is easy to 
show that 



5 K ,K>F aa {K) 



ficlds=0 



xG a (K + K)T (2<1) (K + Ka; Ka; Ka) , (BIO) 



so that 



£<r(*Q 



i I F aa (K)G a (K + K) 

IK 

-(2,1)/ 



x T ( ' ' (K + Ka; Ka; Ka) , (Bll) 



which is the analytic expression for the skeleton graph 
shown in Fig. Ufa). Similarly, we obtain the skeleton ex- 
pansion of the irreducible polarization by differentiating 
Eq. I)B5(I with respect to <fi_R a , 



K ti<P-KcT 5 3K,a8j K+K(T 



ficlds=0 



-i( / G a {K)G a {K + K)T ( ->{K + Ka; Ka; Ka) , 

(B12) 



which is shown diagrammatically in Fig. [5] (b) . Finally, 
applying the operator 



to Eq. 1B5II and sub- 

sequently setting the fields equal to zero we obtain the 
skeleton expansion of the three-legged vertex shown in 
Fig. Etc), 

T {21) (K + Ka;Ka;Ka) = i 
-»C / G a {K')G a {K' + K) 

JK' 

xr (4 ' 0) (if + Ka,K'a;K' + Ka,Ka) . (B13) 

Skeleton expansions for higher-order vertices can be 
obtained analogously from the appropriate functional 
derivatives of Eqs. (jB3)l - (|B7|l . 

APPENDIX C: WARD IDENTITIES 

In this appendix we give a self-contained derivation of 
the Ward identities in Eqs. H4.5fl - I|4.7|) within the frame- 
work of our functional integral approach. Although the 
Ward identity (|4.5|l for the three-legged vertex is well 
known ^M^iS^I it seems that the higher-order Ward 
identities given in Eqs. (|4.6[) and l|4.7|l cannot be found 
anywhere in the literature. Since we are interested in 
deriving infinitely many Ward identities involving the 
vertices r( 2,m ) with two fermion legs and an arbitrary 
number m of boson legs, it is convenient to derive first 
a "master Ward identity" for the generating functional 
for the irreducible vertices, from which we can obtain all 
desired Ward identities for the vertices by taking appro- 
priate functional derivatives. 

Consider the generating functional of the Green's func- 
tion of our mixed Fermi-Bose theory defined in Eq. I|2.18[) . 
which in explicit notation is given by 

G[j,j, J] = — f D$,1>,<p]e- s ®'*'d + G'V + ®*M J ''' / '') . 
J 

(CI) 

If we rewrite the parts of the action involving the 
fcrmionic fields ip and ip in real space and imaginary 
time, the Euclidean action reads [we use again the no- 
tation X = (r, r) introduced in Sec. II V 0] as 

S[tp,ip,(fi] = So[i>, ip] + So[<p] + Si[tp,ip,(p] 
S $,il>] = V [ MX)drMX) 

+ J2j dT \ d ° r Jd D r'i> a {T,r) 

a 

x^(r-r')^(r,r') , (C2) 
SS,1>M = iJ2 I MX)MX)V*{X) , (C3) 

where we have defined the Fourier transform of the dis- 
persion 

dPk 



Ur) 



(2tt 



6c<t e 



ik-r 



(C4) 
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Suppose now that we perform a local gauge transforma- 
tion on the fermion fields, defining new fields tp' and t/r 
via 

MX) = e^ x H' a {X) , i> rT {X) = e- ia ^ x H' a {X), 

(C5) 

where a a (X) is an arbitrary real function. It is easy 
to show that, to linear order in a a (X) 1 the action (|C3|) 
transforms as follows: 

S[e- ia $',e ia ij\<p\ = SW,V, V ] 
+ *W i>'AX)[d T a a (XM a (X) 

a J X 

-i^Jdr Jd D r Jd D r'^(T,r)[a a (T,r)-a a (T,r')} 

x£ CT (r-r')^(r,r') . (C6) 

Using this relation, we see that the invariance of the gen- 
erating functional in Eq. (|Cl|l with respect to the change 
of integration variables defined by Eq. I|C5(I implies, to 
linear order in a a (X\ 



,-S[>P,^,(p] + (j,^)+(i,,j)n+(J,(p) 



= 2o / D bP^>f] 



- J x Mx)[d T a a (x)]Mx) 

+ Y^JdrJd D r Jd D r'MT,r) 

x [a CT (r, r) - a CT (r, r')]£ a (r - r')Vv(T, r') 
+ ti,a*l>)-(#a,j) \ . (C7) 



Taking the functional derivative of this equation with 
respect to a CT (X), this implies in Fourier space, 



s<fig 



sg . sg 

+3k+k<jT^ 3 K <?~ 



K+Ka 



K+Ka 



(C8) 



Expressing this equation in terms of the generating func- 
tional g c = ln(7 of the connected Green's functions and 
the generating functional and T[i(), ip, ip] of the irreducible 
vertices as defined in Eq. (|2.42l) . we obtain 



= 



+1p 



ST 



Ko 



i>K+K. 



8jK*5j K+Ra . 
ST ) 



(C9) 



Alternatively, using the Dyson-Schwinger equation (|B5|I . 



we may rewrite this as 
ST 







5ip R 

(Ck+k,o- 



K 



swg c 



+ I 



K 
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SjKaSj 



K+Ka 



ST 
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K+Ka 



K+Ka 



Ka 



(CIO) 



Equations. i|C9(l and l|C10|) are our "master Ward identi- 
ties" from which we can now obtain Ward identities for 
the vertices by differentiation. For example, taking the 



derivative 



of Eq. (ICTTl we obtain 



iu>Il a (K)-Tl c (T (K) = Q 
where we have defined 



(Cll) 



Tl%(K) = -»C f {^ a -^a)Ga{K)Ga{K + K) 

Jk 

xT (2A) (K + Ka;Ka;Ka) . (C12) 

Equation (jCllf) is a relation between response functions, 
which follows more directly from the equation of continu- 
ity. If we are interested in vertices involving at least one 
fcrmionic momentum and if the momentum transferred 
by the interaction is small, our master Ward identities 
can be further simplified. Then all fermionic momenta 
lie close to a given point kp iCr on the Fermi surface so 
that Eqs. I|C9(1 and IjClOfl become simpler if we assume 
asymptotic velocity conservation. This means that we 
replace under the integral sign 



k+k, a 



k<r 



vf, ct ■ k 



(C13) 



This approximation amounts to the linearization of the 
energy dispersion relative to the point ki? CT on the Fermi 
surface. Using again Eq. (|B5|) . our master Ward identity 
becomes 



= (ioj — VF,a ■ k) 
ST 



ST 



^Ka 



* I i>K+Ko4>Ko 
>K 



+ I 



K 



1> 



Ka 



Sip 



i>K- 



ST 



K 



K+Ka 



Ka 



(C14) 



Differentiating this simplified master Ward identity with 
respect to the fields using the relation (|B8(I as well as 



$(3) r 



5V>KoHKaSlp K+ Ka fio lds=0 

= T {2 ' 1] (K + Ka;Ka;Ka) , 



(C15) 



^k 1 a 5t fK 2 a^KaSlp 



K+Kx+K 2 a 



ficlds=0 



r (2,2) ( K + Kl + k 2(T . Ka . Kl(Ji K 2 a) , (C16) 
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and so on, we obtain the Ward identities for the irre- 
ducible vertices given in Eqs. (|4.51 14.61 14.71) . 

Of course, other Ward identities, e.g., the Ward iden- 
tity for rt 4 ' 1 ) discussed in Ref. El can also be obtained 
from Eq. i|C14() . Note that if the approximation l|C13(l is 
not made, the master Ward identity (|C14p should be re- 



placed by the more general master Ward identity (jClOfl . 
so that the Ward identities l|4.5l 14.61 14.7[> for the vertices 
acquire correction terms. The effect of these correction 
terms on the Ward identities for T^ 2 ' 1 ' and I^ 4 ' 1 ) has very 
recently been studied in a mathematically rigorous way 
by Benfatto and Mastropietro4& 
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